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ABSTRACT 

The  met  lod  of  turbine  performance  prediction  developed  by 
Vavra  and  Eckert  has  been  refined  in  this  analysis  to  realize  more 
of  the  potential  of  the  three-dimensional  calculating  methods. 
Mach  number  and  rotor  tip  clearance  effects  on  blade  outlet  angles 
and  loss  coefficients  have  been  localized  rather  than  averaged  over 
the  blade  height.   An  approximation  for  streamline  curvature  has 
been  used. 

Performance  curves  were  determined  for  two  single  stage  axial- 
flow  turbines  located  at  the  Propulsion  Laboratory  of  the  Naval 
Postgraduate  School.   Test  results  were  available  for  one  of  the 
turbines.   Agreement  between  predicted  and  experimental  performance 
values  was  generally  within  3  per  cent. 


Thesis  by  Robert  G.  Harrison  entitled:   "An  Analysis  of  Single- Stage  Axial- 
Flow  Turbine  Performance  Using  Three-Dimensional  Calculating  Methods" 
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1.   Introduction 

Turbines  form  an  important  part  of  propulsion  systems.   To  opti- 
mize a  design  it  is  necessary  to  know  the  performance  at  off-design 
conditions  as  well  as  the  performance  at  the  design  point.   Since  the 
testing  of  a  prototype  is  very  costly  and  time  consuming,  it  is  of 
great  advantage  to  be  able  to  predict  the  performance  characteristics 
by  means  of  theoretical  methods.   The  more  advanced  a  system  is  the 
more  important  it  becomes  to  improve  the  accuracy  of  these  methods. 
Hence  it  is  necessary  to  base  these  methods  on  the  fundamental  laws  of 
fluid  dynamics  rather  than  on  ru  .e-of -thumb  approximations.   The  latter 
basis  is  possible  only  when  data  on  previous  designs  are  available. 
For  new  and  advanced  configurations,  it  will  be  necessary  to  apply  re- 
fined methods  which  enable  the  designer  and  systems  engineer  to  pre- 
dict the  effect  of  proposed  design  changes. 

The  principal  equations  that  describe  the  flow  properties  in 
turbomachines  are  well  known.   These  same  equations  are  used  as  a  basis 
for  ail  proper  "three-dimensional"  calculating  methods  that  have  been 
cited  in  the  technical  literature.   The  methods  differ  however  in  the 
manner  in  which  these  equations  are  manipulated  and  applied. 

This  thesis  is  concerned  with  the  refinement  of  the  three- 
dimensional  method  of  analysis  developed  by  Vavra  and  Eckert.   Based  on 
the  physical  dimensions  of  particular  test  turbines  that  are  available 
at  the  Turbo-Propulsion  Laboratory  of  the  Naval  Postgraduate  School, 
performance  curves  were  determined  for  these  machines.   Concurrently 
with  the  performance  analysis,  experimental  tests  were  conducted  on 
one  of  these  turbines  so  that  actual  experimental  results  could  be 
used  to  julge  the  accuracy  of  the  proposed  theoretical  performance 
evaluation. 

In  addition  to  his  publications  and  classroom  lectures  which  con- 
stituted the  foundation  for  this  analysis,  Professor  Vavra  was  very 
generous  in  providing  guidance  and  counsel  during  the  period  of  this 
work.   For  this  I  am  greatly  appreciative,   I  would  also  like  to  thank 
Lieutenants  P.  M.  Commons  and  J,  A„  Messegee  for  making  their  experi- 
mental results  available. 
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2.   Basis  for  the  Analysis 

The  two  conservation  equations  that  must  be  satisfied  to  obtain 
a  solution  for  the  flow  in  a  turbine  are  the  equations  of  motion  and 
continuity.  These  equations  are  satisfied  at  stations  between  blade 
rows.  The  method  used  is  that  given  by  Vavra.  Vavra  developed  the 
equations  of  motion  and  continuity  for  absolute  flows.  The  equa- 
tions in  this  form  are  readily  useable  for  the  position  after  the 

stator.   Eckert  later  developed  these  equations  for  relative  flows  to 

2 
be  used  for  rotor  calculations.    Eckert' s  conversion  allows  relative 

flow  quantities  to  be  handled  directly  without  conversion  to  an 

absolute  system.   Eckert' s  approach  also  avoids  iteration  procedures 

to  determine  the  total  enthalpy  after  the  rotor. 

The  assumptions  made  for  the  developement  of  the  equations 

used  in  the  performance  analysis  are: 

1.  An  infinite  number  of  blades  in  each  row  so  that  downstream 
effects  are  not  felt  upstream. 

2.  Axisymmetric  flow  at  the  stations  where  the  equations  of 
motion  are  solved. 

3.  Adiabatic  and  steady  flow  so  that  the  total  enthalpy  along 
any  given  streamline  is  constant  through  the  stator  and  the  relative 
total  enthalpy  is  constant  through  the  rotor. 

4.  All  entropy  changes  are  assumed  to  occur  in  the  blade 
channels  that  are  located  ahead  of  the  stations  where  the  equations 
of  motion  are  satisfied.   Hence  at  the  calculating  stations  the 
flow  is  assumed  to  be  isentropic  along  particular  streamlines. 

With  the  above  assumptions,  the  equations  of  motion  for 
absolute  and  relative  flows,  respecti/ely,  are 

y|-|  =  \7  x  (v*V^   +  TVS 


(1) 

v  HP=  W  x  (7  x  W  +  a  uj)  +  T  vs  (2) 


vavra,  M.  H. ,  Aero-Thermodynamics  and  Flow  in  Turbomachines. 
New  York,  London:   John  Wiley  and  Sons,  Inc.;  1960,  Chapter  16. 
2 
Eckert,  R.  H.,  Performance  Analysis  and  Initial  Tests  of  a 

Transonic  Turbine  Test  Rig  (USNPGS  Thesis,  May  1966),  pp.  149-155. 
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Relative   totei    enthalpy  can  be  written  as 

u       ,         W,2      U,2  _  h         \Nc       Ul  _  L        .  V±tn  _  yl  (3) 

Equivalent  erthalov  i«  oefined  as 

Similar  to  ti..,  the  equivalent  enthalpy  H_  is  constant  along  a 
streamline  for  the  adopted  assumptions.   Equivalent  enthalpy  can 
also  be  written  as  ->     2.   ^ 

HE=h^^  =  K+—  +-ra  (5) 

The  introduction  of  the  equivalent  enthalpy  allows  this  quantity  to 
be  used  for  the  rotor  in  a  manner  analogous  to  the  way  total  enthalpy 
is  used  for  the  stetor. 

For  equations  used  in  this  analy* is,  the  subscripts  refer  to: 

0  -  station  ahead  of  stator 

1  -  station  ahead  of  rotor 

2  -  station  after  rotor 

is  -  isentropic  expansion  from  P 

o 
s  -  isentropic  expansion  from  P 

E 
th  -  the  theoretical  value 

Figure  1  is  a  temperature-entropy  diagram  showing  the  thermo- 
dynamic process  along  a  particular  streamline  for  a  single  stage 
turbine.   In  general,  the  fluid  properties  will  vary  from  stream- 
line to  streamline,   T>-e  methoo  in  which  the  loss  coefficients  are 
applied  is  also  indicated  in  Fig.  1.   The  loss  coefficients  are 
defined  a« 

jfs-  ™ V —  For  the  Stator  (6) 

(.         W*th-Wa  For  the  Rotor  (7) 

The  coordinate  system  that  will  be  used  in  the  analysis  is 
indicates  in  Fig.  2.  This  figure  also  shows  the  general  layout  of 
the  ty|>e  turbine  to  which  the  prediction  performance  analysis  is 
applied. 

ier;  convention  for  the  various,  angles  that  are  needed  in  the 
analysts  {  indicated  in  /i&.   3. 
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The  modification  of  Eq.  2  into  a  form  that  can  be  used  for  the 

analysis  is  given  in  Appendix  A,  Section  1.   The  appropriate  form  of 

Eq.  1  can  be  obtained  from  the  modification  of  Eq.  2  if  the  angular 

velocity  a/   is  set  equal  to  zero.   Other  differences  of  the  final 

equations  derived  from  Eqs.  1  and  2  are  listed  in  Appendix  A, 

Section  1.   Equation  1  can  then  be  written 

)z\  As*l     dot, 


-    ~ —  <~,iN~r<     ■+■       -v — ! —I  —~ — s -  S  M    «,  \ 

x,    -jiNCK'  +   Y,-VA^    dX,     L  Y,av^  'J 

n  2   becomes 


dX,  (8) 


Equation  2  becomes 
dX 

VWY.        W^Y^      %l^   dX^  Lw.iV,2  "   -jclXa 
where: 

y-^L  or  -^        for  Eqs.  8  and  9  respectively 

**■*••     ^        (subscript  m  refers  to  mean  streamline) 

K=5.0  streamline  curvature  factor 

=>r  streamline  displacements  shown  in  Fig.  2 

.  _  L.  +  U 

L-   a  (see  Fig.  2) 

cr 

The  equation  of  continuity  is  used  in  its  non-dimensional  form 
by  introducing  a  flow  function  <£   .   Developement  of   $    is  given 
in  Appendix  A,  Section  2.   In  differential  form  the  equations  for 
absolute  and  relative  flows  can  then  be  expressed  by 


(9) 


4fe&fr-dA}«aA  imtfiir) 


(10) 


and 


fferential  element  of  area  dA  is 

JA-  I  C  oc\r  (12) 
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where: 

2  ■  number  of  blades 

a  =  blade  exit  opening  (see  Fig,  46) 

*  ■  arc.  restriction  coefficient 

Since  £   is  valid  for  isentropic  flow  only,  the  restriction 

factor  *   must  be  introduced  to  correct  the  actual  flow  area  to  an 

effective  area  which  accounts  for  the  restrictions  due  to  the 

boundary  layers  on  both  sides  of  the  flow  channel. 

3 
The  factor  £'   can  be  expressed  by 

In  Eq.  13,  £  *  is  the  boundary  layer  displacement  thickness  and 
h***  is  the  so-called  energy  parameter  defined  as  the  energy  thick- 
ness  divided  by  the  displacement  thickness.   The  tern  — —  repre- 
sents the  loss  that  is  assumed  to  occur  from  the  inlet  to  the 
throat  of  the  blade  channel,  where  ^J   is  the  loss  coefficient  repre- 
senting all  the  losses  across  the  row  of  blades.  The  profile  loss 

4 
coefficient  was  used  by  Eckert   to  represent  the  loss  prior  to  the 

blade  throat.   Percentage  of  the  total  loss  due  to  profile  losses 
will  vary  considerably  depending  on  blade  geometry,  radial  position, 
and  the  incidence  of  the  flow  on  the  leading  edge  of  the  blade. 
Since  secondary  flow  and  tip  clearance  effects  result  in  losses  in 
the  blade  channel,  half  the  total  loss  coefficient  provides  a  better 
average  representation  of  the  losses  in  the  blade  channel  prior  to 
the  throat.   Tie  basis  for  the  developement  of   C   and  H***  as  used 
is  given  in  Appendix  A,  Sections  3  and  4. 

by  multiplying  and  dividing  by  a^  ,„,        Eq,  12  is 

i/>-  :  -      J.X  :   <  .  (14) 


Vavra,  M,  li,,  Problems  of  Fluid  Mechanics  in  Kadial  Turbo- 
machines  (lihode-Saint-Cenese,  Belgium:  Von  Karmin  Institute  for 
Fluid  Dynamics,  1965)  VKI  Course  Note  55b,  pp.  G46-50. 

*n  i 

cckert,  o_p_,  cit .  ,  p.  44. 
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After  integration  Eqs.  10  and  11  become,  respectively, 


vvfr^  /  r 


to 


3 


(15) 


and 


:  w  can  be 
eference 

yvE(T] 


(16) 


The  flowrate  w  can  be  computed  from  the  conditions  ahead  of  the 
stator.   Then  a  reference  flowrate  is  defined  by 


to       >     J      J 

where  W  is    in  square   inches.      Continuity  will  be  satisfied   for 

the  stator  by 

M  -i 

dX         =  wref 

STATOR 


Similarly  for  the  rotor 


¥  4  ^idX 


»w 


rcf 


(17) 
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The  in  luence  of  the  leakage  flow  through  the  radial  tip  clearance 
has  not  been  accounted  for  in  Eq.  19. 

The  element  of  area  between  the  blade  tips  and  the  shroud  is 

Afi=  atrrdr 


The  flow  through  the  tip  clearance  area  is 

-i        Hi 

=  2. 


JX 


""        J   7  ,P 

he  tip  clearance  is 

I        »  w/ 

Jith   th( 


(20) 


(21) 


iP 

""  LEftP^NCE 

Since  the  tip  clearance  is  relatively  small,  the  values  of 

P   ,  T   ,   and  X0  for  the  tip  will  be  used  in  Eq.  21. 
tE   tE        Z 
With  these  assumptions  Eq.  19  can  be  expressed  by 


RcTdK 


vfyret 


Law.*V«]ieeroR 

The  assumptions  used  to  arrive  at  Eq.  22  are  obviously  incorrect  In 


two  respects.   First,  the  flow  represented  by   |    is 


not 
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perpendicular  to  the  tip  clearance  area.   Second,  the  effective  area 
represented  by   c    is  larger  than  that  which  probably  occurs 
because  of  the  relatively  large  boundary  layers  that  exist  on  the 
shroud  and  blade  tips.   The  exact  behavior  of  the  flow  in  the  small 
region  between  the  shroud  and  the  blade  tips  is  impossible  to  predict 
without  further  tests.   However,  it  is  felt  that  the  flowrate  through 
this  space  as  represented  in  Eq.  22  is  too  large  for  the  reasons  just 
mentioned.   Therefore  a  more  accurate  approximation  of  this  flowrate 
will  be  obtained  if  the  last  term  on  the  left  side  of  Eq.  22  is 
divided  by  2,  yielding 


f   Pf. 


ttlHx  W(f(fHX)ripi!ij 


C+»»~c'1r*Ji  ROTOR 

The  tip  clearance  flow  included  in  Eq,  23  can  be  obtained  by  inte- 
gration, 

Introducing  this  expression  into  Eq.  23  gives 


ROTO  P. 


Wr£l 


[aw2^]RaroR  (25) 

3.   Technique  for  Obtaining  Solution 

With  the  equations  of  motion  and  continuity  in  the  forms  given 
by  Eqs.  8,  9,  18  and  25,  a  method  has  been  developed  to  analyze 
single  stage  axial  turbines,  in  particular,  those  available  for  test 
in  the  Turbine  Te6t  Rig  of  the  Turbo-Propulsion  Laboratory  of  the 
Naval  Postgraduate  School,   The  method  of  analysis  predicts  turbine 
performance  for  specified  values  of  inlet  total  pressure,  inlet 
total  temperature,  rotor  speed  and  the  ratio  of  total  inlet  to  static 
discharge  pressure  p  /P2#   This  method  is  similar  to  that  described 
by  Eckert.   However, °Eckert 's  analysis  neglected  some  effects  which 


5Ibid,  Section  3. 
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have  been  accounted  for  in  this  developetnent.   Some  significant 
changes  made  in  this  method  are  listed  below: 

1.  Stator  and  rotor  outlet  angles  for  a  particular  radial 
location  are  computed  using  a  calculated  Mach  number  for  that  loca- 
tion rather  than  an  assumed  Mach  number  or  the  Mach  number  of  the 
mean  streamline. 

2.  Variation  of  loss  coefficients  due  to  changes  in  blade  geo- 
metry in  the  radial  direction  is  accounted  for. 

3.  The  influence  of  rotor  tip  clearance  on  the  rotor  outlet 
angles  and  loss  coefficients  is  concentrated  near  the  tip  of  the 
blade  rather  than  averaging  these  effects  over  the  full  blade  height. 

4.  Fourth  order  polynomials  are  used  to  better  approximate 
the  curves  respresenting  blade  characteristics  as  a  function  of 
radius  and  the  curves  of  rotor  loss  coefficients  as  a  function  of 
incidence, 

5.  Streamline  curvature  effects  have  been  accounted  for  in  the 
solution  of  the  equations  of  motion. 

In  addition  to  the  assumptions  that  were  mentioned  in  Section 
2  for  the  developetnent  of  the  particular  form  of  the  equation  of 
motion,  the  conditions  ahead  of  the  stator  are  assumed  to  be  uniform; 
that  is,  the  total  temperature,  velocity,  and  entropy  are  assumed 
constant  and  the  flow  axial  in  direction.  It  is  realized  that  com- 
pletely uniform  conditions  are  difficult  to  obtain,  but  any  other 
aesumpt  on  would  be  extremely  difficult  to  develop  mathematically. 

Direct  solution  of  the  equations  of  motion  is  not  possible  since 
they  are  nonlinear  in  the  dependent  variable  V,  Likewise,  no  direct 
method  Is  possible  to  satisfy  continuity.  Solutions  of  these  equa- 
tions must  therefore  be  gained  by  making  Initial  assumptions  for  the 
values  of  the  axial  velocities  which  must  be  Improved  by  successive 
Iterations  until  the  equations  are  satisfied.  To  account  for  stream- 
line curvature  and  slope,  a  complete  solution  of  the  flow  through 
the  stator  and  rotor  must  first  be  made  by  neglecting  the  effects 
of  curvature  In  order  to  determine  streamline  locetlons.  Then  the 
iteration  to  account  for  these  effects  may  progress.  These  require- 
ments make  the  use  of  a  highspeed  computer  a  necessity. 
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This  analysis  has  been  programmed  for  the  IBM  360  computer 
using  FORTRAN  IV.  The  program  is  described  in  Appendix  C.   The 
following  paragraphs  set  forth  the  procedural  steps  of  the  program. 
The  equations  are  listed  in  general  form  without  referring  to 
specific  streamline  locations.   In  the  interest  of  clarity,  however, 
some  relationships  will  be  written  in  forms  similar  to  those  used 
in  the  program.   For  example,   ^^Vff0=[T|is(2)/Tf  l*"1    wil1  rePre" 
sent  the  isentropic  relationship  for  the  number  2  streamline. 

Five  streamlines  are  utilized  for  the  analysis  with  the 
number  1  streamline  located  at  the  hub  and  the  number  5  streamline 
located  at  the  tip  as  shown  in  Fig.  2.   The  number  3  streamline  will 
be  used  as  the  mean  streamline,  and  the  radius  of  this  streamline 
will  be  referred  to  as  the  mean  radius.   The  radial  locations  of  the 
streamlines  ahead  of  the  stator  will  be  such  that  the  mass  flowrate 
between  adjacent  streamlines  is  25  per  cent  of  the  total  flowrate. 
Positions  for  the  streamlines  after  the  stator  and  after  the  rotor 
are  initially  assumed.   The  locations  of  streamlines  2,  3,  and  4 
then  vary  during  the  solution  as  necessary  so  that  the  percentage  of 
the  total  flowrate  between  adjacent  streamlines-  does  not  change. 
This  continuity  requirement  will  be  called  streamline  continuity. 

Besides  the  radii,  sufficient  input  information  must  be  used 
to  effectively  reflect  the  physical  characteristics  of  the  stator 
and  rotor  blading.   Some  of  the  physical  properties  are  introduced 
directly;  such  as,  the  number  of  stator  blades,  the  number  of  rotor 
blades,  and  the  rotor  tip  clearance.   The  other  quantities  used  which 
reflect  blade  characteristics  are  throat  opening  dimensions  for  the 
blade  channels,  discharge  angles,  rotor  blade  inlet  angles,  loss 
coefficients,  and  stalling  incidences  for  the  rotor. 

Throat  opening  dimension  "a"  is  a  function  of  radius.  The 
best  method  for  introducing  this  characteristic  into  the  anaylsis 
is  to  enter  the  measured  values  of  "a"  together  with  the  correspond- 
ing radii.  Then,  utilizing  the  method  of  least  squares,  a  fourth 
order  polynominal  curve  is  fitted  through  these  points.  From  the 
resulting  polynomial,  the  value  of  "a"  for  any  radius  required  by 
streamline  continuity  can  be  determined. 
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Discharge  angles  are  nrndicted  by  using  a  combination  of  the 
methods  of  Vavra  and  Ainley.   Outlet  angles  are  first  calculated 
using  the  formula  which  Vavra  established  from  the  experimental  data 
of  Beei .   These  angles  are  then  corrected  for  tip  clearance,  blade 
curvature,  and  Mach  number  effects  with  the  methods  given  by  Ainley. 
Stator  and  rotor  discharge  angles  are  predicted  at  three  radii; 
namely,  the  hub,  mean  radius,  and  tip.   The  method  used  for  calcula- 
ting these  angles  is  explained  in  Appendix  B,  Section  1. 

Values  of  stator  gas  outlet  angles  ex^    for  the  mean  stream- 
line are  determined  tor  Mach  numbers  M.  of  0.5,  0.7,  0.75,  0.8,  and 
1.0.   These  values  are  represented  by  two  parabolic  curves  of  the 
form 

a,=  a+  bM,-  cM*  (26) 

The  first  curve  is  used  for  Mach  numbers  M.  from  0.5  to  0.75  and  the 
second  established  interim  values  of  CX. ,   for  M.  between  0.75  and 
1.0.   From  these  curves  the  flow  angle  a,   for  the  mean  streamline 
cnn  be  found  for  any  Mach  number  M, . 

The  flow  angle  oc,   is  also  a  junction  of  radius  r..  There- 
fore the  changes  of  ex,   ;or  the  hub  and  tip  with  reference  to  the 
mean  radius,  called  A(XH  and  &&<?      respectively,  must  be  used. 
The  flow  angle  at,    <~an  then  be  determined  for  any  r.  by  assuming 
a  linear  variation  ot   OC,    between  the  hub  and  the  mean  radius  and 
between  the  mean  and  the  tip,   With  this  assumption  and  using  the 
Mach  number  M.  in  Eq.  26  corresponding  to  the  number  2  streamline, 
the  flow  angle   or,    for  this  streamline  would  be 

Ainley,  D.  c;,  <md  Mathieson,  G.  C.  R,,  A  Method  of  Performance 
estimation  for  Axial-Flow  Turbines.  Aeronautical  Research  Council, 
I   &  M  No.  2974,  1957.  pp.  3-4. 

Beer,  K,f  Aerodynamic  Design  and  Estimated  Performance  of  a 
Two-Stage  Curtis  Turbine  for  the  Liquid  Oxygen  Turbopump  of  the  M-l 
Engine.  NASA  CR  54764  (AGC  8800-12),  Nov.  19,  1965.  p.  29. 
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The  superscript  **  is  used  with  r-  and  flf   in  Eq.  27  to  indicate 

the  radius  initially  assumed  for  the  mean  streamline  and  the 
computed  for  that  radius.   Equation  27  can  then  be  used  throughout 
the  analysis  even  though  the  radial  location  of  the  mean  streamline 
may  change  due  to  streamline  continuity  requirements.   A  similar 
approach  is  used  to  establish  the  flow  angle  c>     at  the  rotor  dis- 
charge. 

The  rotor  blade  inlet  angles  &£    are  measured  from  the  manu- 
facturing drawings  of  the  blade  profiles.   Using  the  values  ofp^ 
for  the  hub,  mean,  and  tip  streamlines,  a  parabolic  curve  is  deter- 
mined which  gives  J~)         as  a  function  of  X  . 

Loss  coefficients  and  stalling  incidences  are  predicted  by 

8 

using  the  methods  of  Ainley.    For  the  present  method,  the  stalling 

incidence    i   is  defined  as  that  at  which  the  loss  coefficient  is 
s 

twice  the  value  of  the  minimum  loss  coefficient.   Following 

Ainley' s  methods,  loss  coefficients  are  computed  as  a  function  of 

the  ratio  of  flow  incidence  to  stalling  incidence  -Jt-      .      Loss 

L  5 
coefficients  are  also  a  function  of  blade  geometry  or  radius.   Since 

the  stator  has  zero  incidence,  its  loss  coefficient  is  a  function 
of  the  radius  only. 

Stator  and  rotor  loss  coefficients  and  rotor  stalling  incidences 
are  calculated  for  the  hub,  mean,  and  tip  radial  locations.   Potor 
loss  coefficients  £       are  determined  for  values  of   yr    ranging 
from  -2.0  to  1.6.   Curves  of  £g    vs.  ~¥l         are  drawn  for  each  of 
the  three  radial  locations,  and  values  of  &L    ,  along  with  the 
corresponding  nuantities   -7-1   ,  are  used  to  determine  fourth 
order  polynomials  which  approximate  these  curves.   A  similar  proce- 
dur  is  followed  to  obtain  a  fourth  order  polynomial  representing 
stalling  incidence  i   as  a  function  of  radius  f.     .   Sample  calcula- 
tions for  the  prediction  of  stator  and  rotor  loss  coefficients  and 
stalling  incidences  are  contained  in  Appendix  B,  Section  2. 


Ainley,  oj>.  cit.  ,  pp.  4-5. 
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Variation  of  loss  coefficients  with  radial  location  is  accounted 
for  by  assuming  a  linear  variation  of  these  quantities  between  the 
hub  and  the  mean  radius  and  between  the  mean  and  the  tip.   To  demon- 
strate the  procedure  followed  in  computing  rotor  loss  coefficients, 
the  following  example  is  given.   For  a  particular  incidence,  the 
first  step  in  the  determination  of  ^f       for  the  number  2  streamline 
would  be  to  calculate  i  by  using  the  radius  r,  of  that  streamline 
and  the  polynomial  of  the  form  i  ■  i  (r.).  Loss  coefficients  for 
the  hub  and  mean  radius  would  then  be  computed  using  the  resulting 
-4—    in  the  polynomials,  for  these  radii,  of  the  form  %  =  <7C(-^-). 
With  the  assumed  linear  variation  J    (z\      would  be 


%w-%m 


r**-r,(o 


■3;-  x® 


The  reason  for  using  the  superscript  **  on  the  mean  streamline 
value*  Is  the  same  at  previously  mentioned  in  connection  with 
Be.  27. 

For  the  first  approximation, the  Mach  number  M  ahead  of  the 
stator  is  assumed, and  the  static  properties  and  flowrate  at  station 
0  are  round  by 


Te< 

v0 

■frgR-r;     M 
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/: 

(28) 


(>9) 

(30) 

(31) 

(3?) 
(33) 


vV«/4^Vo  (34) 

The  reference  flowrate  W  f  will  ba  used  to  check  continuity  at 
stations  i  and  2. 


2H 


The  next  step  is  to  determine  axial  velocities  after  the 
stator  that  satisfy  the  equation  of  motion.   Total  enthalpy  after 
the  stator  is  constant  by  assumption,  and  streamline  curvature 
effects  are  neglected  at  this  stage  of  the  analysis.   With  these 
conditions,  Eq.  8  simplifies  to 

This  equation  can  be  integrated  to  give 

X 


(36) 


**  /  „t  (37) 

where  X  is  arbitrary  and  //?C  *-s  the  constant  of  integration. 
Using  the  boundary  condition  Y  ■  1.0  at  X  »  1.0,  /t)C  can  be 
found  by 

0=fUt  +  lhCt     on   /aCz=-J~JcJX 


then 


Eauation  38  must  be  expressed  in  a  form  that  can  be  utilized 
in  the  computer.   Expansion  by  infinite  series  yields  ^ 

The  quantities  contained  in  I,  Eq.  36,  must  be  evaluated 
before  proceeding  with  the  solution.   For  the  first  approximation, 

assumptions  are  made  for  the  values  of  Y  ,  M  ,  and  V.  .   The  flow 

m 
angles  (%.     are  then  calculated  by  using  Eqs.  26  and   27.   After  the 

value  of  C[.     has  been  determined  for  each  streamline,  •,^7I     is 

arX  i 

computed.   This  derivative  and  all  others  needed  in  the  analysis  are 
found  by  finite  difference  methods. 
Enthalpy  is  computed  by 

9 
and  the  entropy  term  is  found  by  the  method  of  Vavra. 


9 
Vavra,  M.  H. ,  Aero-Thermodynamics  and  Flow  in  Turbomachines. 

New  York,  London:   John  Wiley  and  Sons,  Inc.,  1960.   pp.  445-447. 
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(38) 


(39) 


(40) 


JL, 


I  - 


Y*  V^ 


(41) 


to  Eq.  28. 

Each  time  a  solution  to  Eq.  36  is  found,  the  new  values  of  Y 
are  then  used  for  the  next  iteration.   After  five  iterations,  0(] 
and   -rr  are  recalculated  using  the  new  stator  exit  Mach  numbers, 
The  Mach  number  for  each  streamline  is  found  by 


V.- 


cosoc. 


'l         'to      agJcp 
With   the  new  values   for     o( 


and 


doc, 


five  more  iterations  are 


made  to  determine  the  corrected  values  of  Y. . 

The  quantities  represented  by  Eqs.  42-44  are  recomputed  after 
satisfying  the  equation  of  motion  and  additional  quantities  deter- 
mined by 


T*%*  \ 


(42) 
(42a) 

(A3) 
(44) 


The  flowrate  through  the  stator  is  computed  next  and  compared 
with  that  required  to  satisfy  continuity.   Reference  flowrate  between 
the  hub  and  each  streamline  is  found  from 

/  / 


(45) 
(46) 
(47) 


(48) 


where 


H»[(i)  -(tP] 


(49) 
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The  "a"  and  "a  "  in  bq.  4b  are  found  from  the  polynomial  which 
represents  throat  opening  as  a  function  of  radius. 

re  .   <*b  is  solved,  the  pressure  ratio  is  compared  with 
the  critical  pressure  ratio,   If  the  critical  pressure  ratio  has 
been  exceeded,  the  flow  is  choked  at  that  radial  location  and  the 
critical  pressure  ratio  is  used  in  calculating    ±>  and    ~ 

The  fraction  of  the  total  flowrate  passing  between  the  hub  and 
each  streamline  is  computed  by 
yy  wi.lX 

Total  flowrate,  as  found  by  the  denominator  of  Eqa  50  multiplied 
hy  -. ..  r  »'..   .  is  then  compared  with  the  reference  flowrate. 
Overall  continuity  is  satisfied  if  the  difference  is  less  than 
0.0002, 

If  required  and  computed  flowrates  are  not  within  this  toler- 
ance, the  axial  velocity  for  the  mean  streamline  is  adjusted  by 

V      -v,     ,  ♦  ;'f^~W;-^j;"i  (51) 

V"  *«.   (N £  .V)   VA I  wl  L  •  O      O.OOQ6? 

is  the  required  reference  flowrate  divided  by  a^.H^V^,  and 

,  is  the  denominator  of  Eq,  50, 
computed 

Solutions  to  the  equation  of  motion  and  continuity  are 
successively  found  until  overall  continuity  is  satisfied, 

The  fractions  of  the  total  flowrate  determined  for  each 
streamline  by  Eq,  50  are  then  compared  with  the  corresponding  flow- 
rate  fractions  ahead  of  the  stator,   Streamline  continuity  is  satis- 
fied if  agreement  is  within  0.002  for  each  streamline.   If  stream* 
line  continuity  has  not  been  satisfied,  the  streamlines  in  error  are 
adjusted  by 


=  XOLD  +fwf     -wf,       i   il 

':  .'/         '^OLO       Tree,.  T cor-  \  J 


JW,  (52> 


equation  52  applies  only  to  stieamlines  2,  3,  and  4  since  W  =  0  at 


the  hub  and  W  ~    1.0  at  the  tip. 
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If  streamline  positions  have  been  adjusted,  new  streamline 
radii  are  found  by 

NEW  ~  X(s)£W  Cu  0LO  (53) 

and  new  values  of  X  for  each  streamline  are  obtained  by 

X=7^  (54) 

^  NEW 

With  the  new  streamline  radii  the  equation  of  motion  and  overall 
continuity  must  be  satisfied  again. 

The  following  values  are  determined  after  streamline  continuity 
is  satisfied: 

U|"(30)(I2)  f>  (55) 

U,=  {f-   U,  (56) 

I'Vu.-U,  (57) 

W 


■1/  Wu, 


I   cos^',  (59) 

TtE=T,  ^2-3-j-Cp  +  EgTE;  <60) 


HE=TttcP  (6i) 

/TtE  N*^"1 

PtE=P)  T/  (62) 


Rotor  blade  inlet  angles  are  determined  for  each  streamline 
by  means  of  the  parabola  which  establishes  pQ      as  a  function  of  X. . 
Incidence  is  found  by 

1  =  ^,-/4  ^ 
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btalling  incidences  and  rotor  loss  coefficients  are  then  deter- 
mined using  the  procdure  described  in  connection  with  Eq.  28.   If 
the  incidence  ratio  -4 —  is  less  than  -2.0  or  greater  than  1.6, 

;       's 

the  value  for  -4—  is  set  equal  to  -2.0  or  1.6,  respectively. 

One  of  the  quantities  necessary  for  the  equation  of  motion  after 

the  rotor  is   -  ba   .  This  quantity  is  separated  into  two  parts, 

d  r^z 

-'L^'9..     represents  the  entropy  gradient  due  to  changes  of  entropy 

across  the  stator  and  referred  to  station  2.     The  term   — — ~ 

d  *a 
represents  the  entropy  gradient  due  to  the  different  entropy  changes 

through  the  rotor.   The  entropy  increase  through  the  rotor  is  com- 
puted by  using  the  corresponding  values  of  the  rotor  in  Eq.  41;  for 
examplejh  would  be  used  instead  of  H. 

Neglecting  streamline  curvature  and  slope,  Eq.  9  can  be  re- 
written 

The  discharge  angle  /^  and  its  derivative  -=p-*  are  found  in  the 
same  manner  as  previously  described  for  the  stator  discharge  angles. 
The  method  for  solving  Eq.  65  is  similar  to  that  used  for  Eq.  36 
with  the  exception  that  iterations  are  carried  out  until  correspond- 
ing values  of  Y  change  by  less  than  0.005yor  un>.il  13  iterations  have 
been  completed.   The  extra  iterations  are  necessary  because  there 
may  be  a  slower  convergence  of  the  values  of  Y  at  station  2.  After 
satisfying  the  equation  of  motion,  values  at  station  2  corresponding 
to  those  represented  by  Eqs.  42-47,  are  computed.   It  should  be 
noted  that  where  absolute  velocity  terms  are  used  in  Eqs.  42-47,  the 
corresponding  equations  for  the  rotor  will  employ  relative  velocities. 
Therefore  M2  is  the  Mach  number  of  the  flow  relative  to  the  rotating 
rotor  blade. 
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Overall  continuity  is  checked  at  the  rotor  discharge  using 
the  same  procedure  utilized  for  the  stator;  however,  the  reference 
flowrate  is  increased  according  to  Eq.  25  to  account  for  the  flow- 
rate  between  the  blade  tips  and  the  surrounding  shroud. 

After  overall  continuity  has  been  satisfied  at  station  2, 
streamline  continuity  is  checked.  In  addition  to  the  functions 
performed  for  the  check  after  the  stator,  there  are  certain  quanti- 
ties that  must  be  adjusted  for  streamline  relocation.  These  are 


^  ■"'/new  \^^2/0Lt.   <:IX2' 


(*fl 


(66) 


(67) 


w  '(^"JL. +  TtJ2,  Cx— *  M  (68) 

A  solution  exists  when  the  equations  of  motion,and  overall 
and  streamline  continuity  have  been  satisfied;  however,  this  solu- 
tion has  neglected  streamline  curvature  and  streamline  slope.  To 
obtain  a  solution  which  accounts  for  streamline  curvature  effects, 
the  terms  in  Eqs.  8  and  9  that  contain  %r    or  LR     must  be  in- 
cluded when  solving  the  equations  of  motion. 

The  streamline  displacements   >'   and  DR,  shown  in  Fig.  2, 
are  computed  for  each  streamline  by 

Sr=r,  -  k£^  (70) 

and 

AR-r.-r,  (7i) 

Equation  70  Is  based  on  the  assumption  that  the  radii  of  any  stream- 
lint;  at  stations  0  and  2  are  not  greatly  different.  Therefore  the 
cosine  of  the  angle  between  the  line  connecting  these  points  and  the 
axis  of  the  machine  is  approximately  equal  to  unity.  The  length  L, 
also  shown  In  Fig.  2,  is  taken  to  be  half  the  distance  from  0.1  inch 
ahead  of  the  stator  to  0.1  inch  after  the  rotor. 
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After  a  solution  is  obtained  which  accounts  for  streamline 

P 
curvature  effects,  the  resultant  pressure  ratio  —=2 —  is 

P  a    P* 

compared  with  the  -£—    specified  initially.   Ideally    to   for 

each  streamline  will  be  the  same.   However,  computer  solutions  for 

this  quantity  may  vary  slightly  from  streamline  to  streamline. 

Therefore  a  mass-flow-weighted  value  of  ~£ —  ,  called  (—^ — )  •  , 

is  found  by 


(72) 


Z    3- 

C-l 

If  the  specified  — ^-  and  (——-). differ  by  more  than  0.0003,  the 
Mach  number  of  the  flow  ahead  of  the  stator  is  properly  adjusted 

and  another  solution  is  found. 

P 
After  Che  iterations  for  — i^—  have  been  completed,  additional 

P& 
quantities  are  determined  for  each  streamline, using 

AH  =  H,-H£--(u1Vo,-U,Voi)gy  (73) 

Overall  efficiency  is  then  computed  by 

f"  Alu.-Tfc-T.fr  (76) 

The  ideal  change  in  enthalpy   Ahi^  is  the  isentropic  enthalpy 
drop  from  the  total  inlet  pressure   F\0   to  the  static  discharge 
pressure  F  .   Equation  76  is  used  to  compute  the  efficiency  of  a 
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single  stage  turbine  because  the  kinetic  energy  leaving  the  rotor 

cannot  be  utilized.   The  efficiency  defined  by  En.  76  will  be 

referred  to  as  total-static  efficiency. 

Theoretical  degree  of  reaction  r*  and  head  coefficient  k.   are 

is 

given  by 


r  = 


and 


For  turbine  performance  curves  it  is  desirable  to  obtain  the 

mass-flow-weighted  values  of  efficiency  (     ).,  head  coefficient 

w 

(k.  ).,  theoretical  degree  of  reaction  (r*).,  horsepower  (HP).,  and 
is'w'  &  v  'w'  v  v   7w' 


<V*- 


The  last  two  nuantities  are  found  from 


(alILIh 


(77) 


(78) 


(80) 


'"^'W   "        SSO  (79) 

and 

M*  -  ^¥^ 

The  mass-flow-weightedZlA/ ,    called    (£   H) .  ,    as  well   as   (  7?     )., 

(k.  ).,  and  (r*).  are  computed  using  enuations  similar  to  Eo.  72. 
is  w         w        r 

Referred  values  are  obtained,  following  NASA  practice. 

For     £       -    1.4: 

///*  =      (HP)*t  (81) 

f       fe   S 

M<>rer    =    (J^U^-  (82) 

o 

WM+f  =  -4^-  (83) 

VVVe/>    =     W  SW  (84) 

x/  V  (84a) 
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where: 

g-  fto     _     Tto_ 
Tsto.  ~  3/8.4 

°  ~  Psro.    "   M.7 

(85) 

(86) 

A.   MOD  I  and  MOD  II  Turbines 

The  method  of  analysis  as  presented  was  used  to  determine 
performance  curves  for  the  so-called  MOD  I  and  MOD  II  turbines. 
Both  turbines  are  single  stage  axial-flow  machines.   Experimental 
tests  were  conducted  on  the  MOD  II  turbine  by  Commons  and  Messegee 
and  are  described  in  Fefs.  4  and  6.   The  test  results  are  plotted 
with,  appropriate  predicted  performance  curves. 

The  so-called  MOD  I  turbine  was  designed  for  free-vortex  flow 
and  has  highly  twisted  blades.   Outer  diameter  of  this  turbine  is 
9.898  inches.   The  hub  diameters  of  the  stator  inlet  and  rotor  dis- 
charge are  6.930  and  5.970  inches,  respectively.   The  stator  con- 
tains 13  blades  and  the  rotor  22  blades.   Blades  of  the  MOD  I  tur- 
bine are  generally  thin.   The  stator  and  rotor  profiles  used  to 
predict  outlet  angles  and  loss  coefficients  are  shown  in  Fig.  4. 

The  MOD  II  turbine  is  approximately  the  same  size  as  the  MOD  I, 
but  its  blading  is  considerably  different.   The  blades  of  the 
MOD  II  turbine  are  thick  with  blunt  leading  edges  and  constant 
profiles  over  the  blade  height.   Outer  diameters  of  the  MOD  II 
turbine  stator  and  rotor  are  9.701  and  9.836  inches,  respectively. 
The  stator  has  a  hub  diameter  of  6.796  inches,  and  the  hub  diameter 
of  the  rotor  is  6.598  inches.   There  are  19  stator  blades  and  18 
rotor  blades.   Stator  and  rotor  blade  profiles  for  the  MOD  II  tur- 
bine are  shown  in  Fig.  5. 

Throughout  the  remainder  of  this  thesis  the  MOD  I  and  MOD  II 
turbines  will  be  referred  to  simply  as  MOD  I  and  MOD  II. 

The  minimum  throat  opening  "a"  of  the  blade  channels  is  a  very 
critical  dimension.   Slight  variations  of  this  quantity  have  a  con- 
siderable effect  on  turbine  performance.   Since  this  quantity  is  so 
sensitive,  values  of  "a"  measured  from  the  actual  hardware  were  used 
for  the  analysis  rather  than  those  obtained  from  the  manufacturing 
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drawings.   Then  errors  due  to  manufacturing  will  not  be  a  factor  in 
comparison  of  predicted  and  experimental  results.   Figure  6  shows 
the  throat  openings  "a"  as  a  function  of  radius  for  the  stator  and 
rotor  blading  of  both  turbines. 

Predicted  stator  outlet  angle  0(^   as  a  function  of  radius  r- 
for  the  MOD  I  is  plotted  in  Fig.  7.   The  assumed  linear  variation  of 
X .       between  the  hub  and  the  mean  radius  and  between  the  mean 
and  the  tip  is  readily  apparent  in  this  figure.   Figure  8  shows  the 
predicted  variations  of   \(  with  Mach  number  M1  for  the  MOD  I. 
Figures  9  and  10  are  the  corresponding  plots  for  the  MOD  II. 

Relative  discharge  angles  j8z   were  computed  for  two  radial  tip 
clearances  k  for  each  turbine;  namely,  0.020  and  0.033  inches  for 
the  MOD  I,  and  0.015  and  0.033  inches  for  the  MOD  II.   The  predicted 
flow  angles  z.      are  plotted  in  the  same  manner  as  previously 
described  for  the  flow  angles   c\  (  .   Figures  11  and  12  show  /£■_ 
as  a  function  of  radius  r-  and  as  a  function  of  Mach  number  M-, 
respectively,  for  the  MOD  I,  where  M_  refers  to  the  Mach  number 
of  the  flow  relative  to  the  rotor.   Figures  13  and  14  show  the 
corresponding  plots  for  the  MOD  II.   The  predicted  effect  of  radial 
tip  clearance  on  the  discharge  angles  pz     can  be  seen  in  Figs.  11  - 
14. 

Stator  loss  coefficients  li1        were  computed  for  the  radial 
locations  corresponding  to  the  hub,  mean  radius,  and  tip.   Figure 
15  shows  'ff         as  a  function  of  radius  r   for  both  the  MOD  I  and 
MOD  II.   The  straight  lines  in  this  figure  between  the  values  of 

at  the  hub  and  mean  radius^and  between  the  mean  and  the  tip, 
reflect  the  assumed  linear  variation  of   ^c   in  these  regions. 

Variation  of  rotor  blade  inlet  angle  A.,     with  radius  r  is 
plotted  for  both  turbines  in  Fig.  16.  The  difference  between  the 
untwisted  and  the  free-vortex  blades  is  easily  noted  in  this  plot. 
Also  shown  in  this  figure  are  curves  representing  the  variation  of 
stalling  incidence  i  with  r..  .   The  change  of  the  MOD  I  rotor 
blade  profiles  with  radius  is  reflected  by  a  considerable  variation 
of  i  whereas  just  the  opposite  is  true  for  the  MOD  II. 
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Curves  for  the  MOD  I  showing  predicted  rotor  loss  coefficients 
Uf  as  a  function  of  incidence  ratio  -4 —  for  the  hub,  mean  radius, 

and  tip  are  shown  in  Fig.  17.   Since  the  loss  coefficients  between  th( 
mean  radius  and  the  tip  are  dependent  on  tip  clearance,  there  are  two 
curves  for  the  tip.   One  curve  holds  for  the  tip  clearance  of  0.020 
inches;  the  other  is  for  the  larger  tip  clearance  of  0.033  inches. 
For  negative  incidence  ratios  between  -0.5  and  -2.0  the  curves  for 
the  tip  are  estimations.   This  was  necessary  because  the  computations 
by  the  method  shown  in  Appendix  B,  Section  2,  gave  unrealistically 
low  values  of  ^L  as  the  flow  inlet  angle  A      approached  -90  . 
The  situation  is  more  easily  understood  when  it  is  noted  that  the 
blade  inlet  angle  is  -33.7  and  the  stalling  incidence  is  37  at 
the  rotor  blade  tip.   Figure  17  shows  that  the  loss  coefficients  for 
the  hub  (r=3.300  in.)  are  relatively  large.   The  loss  coefficients 
at  the  hub  are  larger  for  most  incidence  ratios  than  those  at  the 
tip  for  a  tip  clearance  of  0.020  inches.   The  larger  value  of 
at  the  hub  reflects  the  higher  losses  that  are  associated  with  an 
impulse  type  blade.   It  can  be  seen  in  Fig.  4  that  the  blade  shape 
varies  from  a  reaction  type  profile  at  the  tip  to  an  impulse  type 
profile  at  the  hub. 

Loss  coefficients  for  the  MOD  II  rotor  are  plotted  in  Fig.  18 
for  tip  clearances  of  0.015  and  0.033  inches.   The  predicted 
similarity  of  the  curves  for  the  hub  and  mean  radius  is  to  be  ex- 
pected since  the  blading  differs  only  in  solidity.   Although  the 
blade  profile  is  the  same  at  all  radii,  the  losses  due  to  tip 
clearance  result  in  larger  loss  coefficients  for  the  tip  profile. 

For  convenience,  the  blade  properties  used  for  calculating  the 
MOD  II  rotor  loss  coefficients  were  those  at  the  hub,  mean  and  tip 
radii  of  the  rotor  discharge.   Since  the  annulus  area  at  the  rotor 
discharge  is  larger  than  the  annulus  area  at  the  stator  discharge 
and  since  the  flow  incidence  is  a  significant  parameter  for  the 
rotor  loss  coefficients,  it  would  have  been  more  appropriate  to  use 
the  blade  characteristics  at  the  hub,  mean,  and  tip  radii  of  the 
rotor  inlet.   However,  the  error  is  insignificant  because  the  blade 
profile  does  not  change  along  the  radius. 
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It  may  be  noted  from  Fig.  4  that  the  minimum  radius, for  which 
a  blade  profile  is  given, is  3.597  inches  whereas  the  radius  at 
the  hub  of  the  MOD  I  stator  ex^t  is  3.300  inches.   The  outlet  angle 
and  loss  coefficient  for  the  hub  were  found  by  extrapolation, using 
the  values  computed  for  the  radii  of  4.125  and  3.597  inches  and 
assuming  a  linear  variation  of  these  quantities  with  radius.   The 
relative  outlet  angle  for  the  hub  of  the  MOD  I  rotor  was  found  in 
a  similar  manner. 

Performance  curves  for  the  two  turbines  were  determined  for 
the  rotor  tip  clearances  mentioned  earlier.   The  total  inle  t  to 
static  discharge  pressure  ratios  investigated  were  1.30,  1.40, 
1.50,  and  1.60,  with  the  exception  that  pressure  ratios  of  1.31  and 
1.51  were  used  for  the  MOD  II  with  0.033  inch  rotor  tip  clearance. 
These  pressure  ratios  agree  more  closely  with  those  experimentally 
investigated  by  Commons  and  Messegee. 

The  exial  distances  L  used  for  the  determination  of  the 
curvatures  depend  on  the  axial  clearance  between  the  stator  and 
rotor  as  well  as  on  the  blade  geometries.   Axial  clearances  of  0.4 
and  1.0  inches  were  used  for  the  analyses  of  the  MOD  I  and  MOD  II, 
respectively, 

Curves  representing  the  performance  of  the  MOD  I  are  plotted 
in  Figs.  19  through  26.   Performance  values  plotted  are  mass-flow- 
weighted  values  unless  stated  otherwise.   Figure  19  shows  referred 
flowrate  as  a  function  of  referred  RPM.   The  increase  in  flowrate 
due  to  an  increase  in  rotor  tip  clearance  can  be  seen  in  this 
figure.  Although  the  flowrate  is  greater  for  the  larger  tip 
clearance, the  toi que  developed  is  greater  at  the  smaller  tip 
clearance.   The  decrease  in  torque  for  the  larger  tip  clearance 
results  from  the  increased  losses  and  decrease  in  turning  angle  of 
the  flow  through  the  rotor  near  the  tip.   The  predicted  effect  of 
tip  clearance  on  torque  is  shown  in  Fig,  20  where  referred  moment 
is  plotted  versus  referred  RPM. 

The  variation  of  total-static  efficiency  with  referred  RPM  for 
the  two  tip  clearances  can  be  seen  in  Figs.  21  and  22.   The  referred 
RPM  at  which  maximum  efficiency  occurs  increases  when  the  total 
inlet  to  static  discharge  pressure  ratio  is  increased.   Blade  losses 
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and  the  kinetic  energy  of  the  flow  leaving  the  rotor  affect  the  total- 
static  efficiency.   As  pressure  ratio  is  increased  the  absolute  veloc- 
ity leaving  the  stator  and  the  relative  velocity  leaving  the  rotor 
increase.   Therefore  the  peripheral  speed  of  the  rotor  must  increase 
to  obtain  conditions  where  the  absolute  velocity  leaving  the  rotor 
is  in  an  axial  direction  and  where  the  relative  flow  ahead  of  the 
rotor  has  zero  incidence.   The  RPM  where  the  flow  has  zero  incidence 
on  the  rotor  will  not  necessarily  be  that  at  which  the  absolute 
velocity  leaving  the  rotor  is  axial.  At  any  RPM, flow  incidence  and 
absolute  discharge  angles  vary  from  streamline  to  streamline, and  the 
above  statements  refer  to  mass -flow-weigh ted  values. 

It  may  be  noted  in  Figs.  21  and  22  that  the  peak  total-static 
efficiency  decreases  somewhat  as  the  pressure  ratio  increases.   At 
higher  pressure  ratios  the  ratio  of  kinetic  energy  leaving  the  rotor 
to  the  work  done  on  the  rotor  increases.   Since  the  kinetic  energy 
leaving  the  rotor  is  lost  energy  for  a  single  stage  turbine,  the 
total-static  efficiency  declines.   The  effects  of  pressure  ratio  and 
tip  clearance  on  efficiency  can  be  seen  in  Fig.  23  where  total-static 
efficiency  is  plotted  as  a  function  of  the  isentropic  head  coefficient. 

The  variation  of  referred  power  with  referred  RPM  can  be  seen  in 
Fig.  24.   Peak  power  does  not  occur  at  the  same  referred  RPM  at 
which  peak  efficiency  occurs.   The  peak  referred  power  occurs  at  the 
referred  RPM  where  the  product  of  total-static  efficiency  and  referred 
flowrate  is  greatest. 

Theoretical  degree  of  reaction  is  plotted  as  a  function  of  isen- 
tropic head  coefficient  in  Fig.  25.   It  may  be  noted  that  the  theo- 
retical degree  of  reaction  increases  with  increasing  pressure  ratio  and 
decreases  with  increasing  tip  clearance  for  any  given  isentropic  head 
coefficient.   The  predicted  effects  are  considerably  different  from 
the  results  of  the  radial  turbine  tests  conducted  by  Riley. 
Riley  found  that  theoretical  degree  of  reaction  was  independent  of 
pressure  ratio  and  axial  clearance  for  radial  turbines. 


10Riley,  M.  W.,  The  Effect  of  Axial  Clearance  on  the 
Performance  of  a  Dual  Discharge  Radial  Turbine  (USNPG  Thesis, 
December  1966^ ,  p.  70. 
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Performance  values  for  each  streamline  are  obtained  from  the 
computer  solution.   However,  plots  utilizing  values  for  each  stream- 
line would  be  difficult  to  analyze.   The  deviation  of  the  hub  and 
tip  values  from  that  of  the  mass- flow-weighted  average  may  be  seen 
in  Flo.  26.   In  this  figure, hub,  tip,  and  mass-flow-weighted  values 
of  theoretical  degree  of  reaction  are  plotted  as  functions  of 
referred  RPM  for  a  total  inlet  to  static  discharge  pressure  ratio  of 
1.40. 

Performance  curves  for  the  MOD  II  corresponding  to  those  pre- 
sented for  the  MOD  I  are  plotted  in  Figs.  27  through  39.   Additional 
plots  have  been  used  for  the  MOD  II  because  of  the  inclusion  of 
experimental  results.   An  axial  clearance  of  1.0  inches  was  used  for 
the  theoretical  prediction.   Therefore, only  experimental  data  for 
that  axial  clearance  are  shown. 

Comments  made  concerning  the  performance  curves  of  the  MOD  I 
are  applicable  to  the  MOD  II  performance  curves  also.   Differences 
in  the  performance  of  the  two  turbines  will  be  discussed  later. 

Plots  of  the  variation  of  referred  flowrate  with  referred  RPM 
are  shown  in  Figs.  27  and  28  for  two  rotor  tip  clearances.   The 
quantitative  values  as  well  as  the  curve  shapes  agree  well  with  the 
experimental  data.   The  maximum  difference  between  predicted  and 
experimental  referred  flowrates  occurs  at  a  pressure  ratio  of  1.51 
for  a  tip  clearance  of  0.033  inches.   There  are  two  experimental 
points  for  this  pressure  ratio  and  tip  clearance  that  differ  from 
the  predicted  curve  by  about  2  per  cent.   The  predicted  and  experi- 
mental values  for  all  pressure  ratios  and  tip  clearances  have  an 
average  difference  of  less  than  1  per  cent. 

Figures  29  and  30  show  curves  of  referred  moment  versus 
referred  RPM.   The  trends  expressed  by  the  predicted  curves  ere  in 
excellent  agreement  with  the  experimental  data.   Although  the  quanti- 
tative agreement  between  theoretical  and  experimental  values  is  very 
good  for  three  of  the  curves,  the  experimental  torque  is  generally 
lower  than  the  predicted  torque.   The  average  difference  between  pre- 
dicted and  experimental  values  is  about  3  per  cent. 
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Figures  31  through  34  show  total-static  efficiency  as  a  func- 
tion of  referred  RPM.   The  shapes  of  the  predicted  curves  generally 
agree  well  with  the  experimental  results.   However,  there  is  an  indica- 
tion that  experimental  efficiencies  decrease  more  rapidly  at  high  RPM 
than  is  predicted  by  the  theoretical  curves.   In  the  high  RPM  region 
the  upper  part  of  the  rotor  blade  has  a  large  negative  flow  incidence. 
In  the  prediction  analysis  when  the  incidence  ratio     had  a  value 
less  than  -2.0,  the  value  of  -2.0  was  used  for  computing  loss  coeffi- 
cients.  This  limitation  may  be  the  reason  that  the  predicted  effi- 
ciencies in  the  high  RPM  range  do  not  decrease  as  rapidly  as  the  test 
data  indicate. 

The  quantitative  agreement  between  predicted  and  experimental 
efficiencies  varies  considerably  between  different  pressure  ratios 
and  tip  clearances.   There  are  two  data  points  at  a  pressure  ratio 
of  1.31  and  tip  clearance  of  0.033  inches  where  the  experimental 
efficiencies  are  over  five  points  below  the  predicted  values.   At  a 
pressure  ratio  of  1.50  and  a  tip  clearance  of  0.015  inches  the 
average  difference  between  predicted  and  experimental  efficiencies 
is  1.5  points.   Giving  equal  weight  to  all  experimental  values  the 
average  difference  between  experimental  and  predicted  values  is 
2.6  points. 

The  calculations  gave  a  decrease  in  efficiency  by  about  two 
points  as  the  tip  clearance  was  increased  from  0.015  to  0.033 
inches.   The  decrease  in  experimental  efficiencies  for  the  increased 
tip  clearance  varied  with  the  different  pressure  ratios.   However, 
the  average  decrease  in  efficiency  is  close  to  the  predicted  decrease. 

Total-static  efficiency  as  a  function  of  isentropic  head  co- 
efficient k.   is  shown  in  Fig.  35  for  pressure  ratios  of  1.40  and 
1.60.   Figure  38  shows  degree  of  reaction  r*  versus  k   .   No  experi- 
mental data  are  plotted  in  these  figures  because  experimental  values 
of  mass-flow-weighted  r*  and  k   were  not  available. 

Plots  of  referred  power  as  a  function  of  referred  RPM  are  shown 
in  Figs.  36  and  37.   The  comments  made  earlier  concerning  the  referred 
moment  plots  apply  to  these  curves  also, since  the  turbine  power  is 
proportional  to  the  product  of  torque  and  RPM. 
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Hub,  tip,  and  mass-flow-weighted  values  of  theoretical  degree 
of  reaction  are  plotted  in  Fig.  39  for  a  pressure  ratio  of  1.40. 
The  experimental  values  for  the  hub  and  tip  er^   also  plotted  in  this 
figure.   The  trend  of  the  experimental  data  is  in  close  agreement 
with  the  predicted  trend;  however,  the  quantitative  agreement  is  poor, 
especially  at  the  hub.   The  decrease  in  degree  of  reaction, as  the  tip 
clearance  is  increased ; is  considerably  greater  than  predicted.   Both 
predicted  and  experimental  results  indicate  that  an  axial  flow  tur- 
bine diffeis  from  a  radial  turbine  inasmuch  as  the  theoretical  degree 
of  reaction  changes  if  rotor  tip  clearance  is  changed. 

Figure  40  shows  the  absolute  flow  angles  CX.,  and   0(  j    as  a 
function  of  radius  for  the  stator  and  rotor  outlets,  respectively, 
of  the  MOD  II.   Experimental  data  were  obtained  for  the  same  opera- 
ting conditions;  namely,  P  /P  =1 .40,k=0. 015  inches,  RPM//^  =13,394. 

co  L 
The  experimental  and  predicted  stator  outlet  angles  have  an  average 

difference  of  less  than  half  a  degree  over  the  blade  height.   Pre- 
dicted and  experimental  agreement  for  the  absolute  flow  angles  at  the 
rotor  outlet  is  not  as  close.   If  the  experimental  point  near  the  hub 
is  neglected,  the  average  difference  between  predicted  and  experimen- 
tal values  of  ">  a  is  about  3.5  degrees.   The  experimental  flow 
angles   A    are  larger  than  the  predicted  angles  over  most  of  the 
blade  height,  which  explains  why  the  predicted  torque  was  generally 
larger  than  the  measured  torque. 

Velocity  distributions  at  the  stator  and  rotor  exits  of  the 
MOD  II  are  shown  in  Fig,  41  for  the  operating  conditions  listed  in 
the  preceding  paragraph.   The  experimentally  found  stator  discharge 
velocities  varied  only  slightly  more  from  hub  to  tip  than  predicted. 
The  average  difference  between  predicted  and  experimental  stator  exit 
velocities  was  about  1  1/2  per  cent. 

The  experimentally  determined  absolute  rotor  exit  velocities, 
which  are  shown  in  Fig.  41,  are  considerably  larger  than  the  pre- 
dicted velocities  except  at  the  hub  and  tip.   Although  experimental 
points  could  be  taken  only  at  distances  greater  than  0.16  inches  from 
the  inner  diameter  of  the  shroud,  the  trend  of  the  data  points  indicates 
that  the  experimental  velocities  are  less  than  predicted  at  the  tip. 
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Near  the  hub  the  experimental  velocities  were  less  than  predicted. 
The  decrease  in  experimental  velocities  near  the  hub  and  tip  is 
the  result  of  separated  flow  in  these  regions.   The  reader  is 
referred  to  Kef.  6  for  additional  information  concerning  the  experi- 
mental data  and  for  photographs  of  the  rotor  showing  indications  of 
separation. 

Experimental  relative  rotor  discharge  angles  £?    were  deter- 
mined for  the  known  rotor  speed  from  the  measured  values  of  absolute 
rotor  outlet  angles  and  velocities,  which  are  plotted  in  Figs.  40 
and  41,  respectively.   The  experimental  and  predicted  values  of  /^ 
are  shown  in  fig.  42.   The  magnitudes  of  the  experimental  values 
of   y.  near  the  hub  are  larger  than  predicted.   Over  approximately 
the  outer  three  fourths  of  the  rotor  blade  height  the  magnitudes  of 
the  experimental  values  of  £•       are  less  than  predicted.   At  a 
radius  6f  4.4  inches,  for  instance,  the  experimental  value  of  & 
is  about  -57  degrees  whereas  the  predicted  value  is  about  -67 
degrees.   The  average  difference  between  experimental  and  predicted 
values  of  ~f.      is  about  6.5  degrees.   The  lower  values  of   •r-'/   pre- 
dicted near  the  hub, provide  some  compensation  in  the  overall  performance 
for  the  high  values  of  ,:'       predicted  for  the  outer  part  of  the  blade. 

The  MOE  II  was  not  designed  to  achieve  the  highest  possible 
efficiency.   The  objectives  of  the  design  were  to  investigate  the 
effects  of  b  unt  untwisted  bladings  of  constant  profile.   This  type 
of  blading  has  the  following  advantages: 

1.  Blade  cooling  passages  are  easily  accommodated. 

2.  A  wider  operating  RPM  range  is  posible  for  a  specified  tur- 
bine efficiency  variation, 

3.  Constant  profile  blades  can  be  manufactured  more  economi- 
cally than  twisted  blades,  especially  if  exotic  high  strength  materials 
are  needed  for  elevated  temperatures. 

The  following  discussion  concerning  comparison  of  the  MOD  I  and 
II  is  based  entirely  on  the  predicted  performance  of  these  tur- 
bines using  the  methods  described  in  this  thesis.   The  reader  should 
be  aware  that  Ainley's  methods  for  predicting  loss  coefficients  were 
not  developed  for  use  with  rotor  blades  having  blunt  leading  edges. 
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Therefore  the  accuracy  of  the  rotor  loss  coefficients  predicted 

for  the  MOD  II  is  probably  not  as  good  as  for  the  MOD  I. 

Figure  43  shows  the  predicted  axial  velocity  ratios  as  a 

function  of  radius  ratio  at  the  stator  exit  and  at  the  rotor  exit 

for  the  MOD  I  and  the  MOD  II.   The  curves  are  for  the  conditions 

that  occur  at  the  maximum  efficiency  of  each  turbine.   The  axial 

velocities  of  the  MOD  I  are  nearly  constant  over  the  blade  height  at 

the  stator  exit  and  at  the  rotor  exit.   The  blading  of  this  turbine 

was  originally  designed  for  free-vortex  flow  by  assuring  uniform 

axial  velocity  components  from  hub  to  tip.   From  Fig.  43  it  can  be 

seen  that  this  condition  is  only  approximately  satisfied  since  loss 

variations  and  curvature  influences  produce  slight  variations  fi 

the  assumed  distribution.   Axial  velocities  for  the  MOD  II  decrease 

from  the  hub  to  the  tip  at  the  stator  exit, and  increase  from  the  hub 

to  the  tip  at  the  rotor  exit, since  the  bladings  of  this  turbine  have 

constant  profiles.   The  variation  of  axial  velocities  for  the  MOD  II 

is  more  pronounced  at  the  rotor  exit  than  at  the  stator  exit.   At 

vVA: 
the  rotor  exit   — varies  from  0.735  at  the  hub  to  1.195  at  the 

tip. 

The  maximum  predicted  total-static  efficiency  of  the  MOD  I  at 
a  tip  clearance  of  0.033  inches  is  about  84  per  cent.   The  MOD  II 
has  a  maximum  predicted  efficiency  at  that  tip  clearance  of  about  80 
per  cent.   The  difference  in  predicted  peak  efficiencies  is  cvo   to 
different  factors.   The  stator  loss  coefficients  are  hi}:he   for  the 
MOD  I  than  for  the  MOD  II,  whereas  the  loss  coefficients  for  the 
MOD  II  rotor  at  zero  incidence  are  tomtwhat  larger  than  those  for 
the  MOD  I  at  zero  incidence.   An  accurate  comparison  of  loss  co- 
efficients is  more  involved  for  the  rotor  than  for  the  stator. 
Because  of  its   tw'Jsted  blades  ?the  MOD  I  rotor  has  nearly  zero  inci- 
dence at  all  points  along  the  blade  height  in  the  vicinity  of  the 
optimum  efficiency.   Only  at  one  radius  of  the  MOD  II  rotor  blade, 
however,  is  the  incidence  zero.   The  incidence  angles  at  larger 
radii  are  negative  and  at  smaller  radii  they  are  positive.   There- 
fore part  of  the  blade  is  always  operating  at  a  loss  coefficient 
larger  than  the  minimum.   Of  probably  even  greater  significance 
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for  the  efficiency  decrease  is  the  kinetic  energy  of  the  gas  leaving 
the  turbine.  Minimum  kinetic  energy  is  lost  when  the  flow  is  axial 
in  direction.   The  design  of  the  MOD  I  is  such  that  at  all  points 
along  the  blade  height  the  absolute  velocity  at  the  rotor  exit  is 
nearly  axial  in  the  vicinity  of  the  optimum  efficiency.   The  MOD  II 
has  only  a  small  radial  portion  of  the  rotor  blade  where  the  dis- 
charge angle  (X^  is  zero.  At  radii  greater  than  the  radius  where 
0(.  2l   is  zero,  the  absolute  flow  discharge  angle  is  positive  and  at 
smaller  radii  it  is  negative.   Therefore  the  kinetic  energy  leaving 
the  rotor  of  the  MOD  H  is  ,reater  than  that  of  the  MOD  I  when  both 
turbines  are  operating  at  peak  efficiency. 

A  comparison  of  the  total-static  efficiency  versus  isentropic 
head  coefficient  curves  for  the  MOD  I  (Fig.  23)  and  MOD  II  (Fig.  35), 
shows  a  larger  variation  of  efficiency  with  head  coefficient  for  the 
MOD  I.   For  example,  at  k=0.033  inches  and-—- =1.60,  the  MOD  I 
efficiency  decreases  about  11.5  points  as  the  head  coefficient  in- 
creases from  3  to  7  whereas  the  MOD  II  efficiency  decreases  only  about 
7  points  for  the  same  change  in  head  coefficient.   This  difference  in 
efficiency  variation  results  from  the  effects  of  blade  twist  mentioned 
earlier. 

To  investigate  the  influence  of  the  streamline  curvatures  on  the 
predicted  performance  of  the  MOD  I  and  MOD  II  the  presented  flow 
equations  were  solved  for  an  axial  length  of  the  bladings  L  of  9  x  10 
inches.  Increasing  L  to  this  value  has  the  same  effect  as  neglecting 
streamline  curvature  effects  as  can  be  seen  from  Eqs.  124  and  125. 
Differences  between  performance  values  neglecting  streamline  curvature, 
and  those  where  curvature  effects  were  accounted  for.were  less  than 

0.2  per  cent  for  both  turbines. 

5.   Conclusions  and  Recommendations 

Effects  of  streamline  curvature  were  found  to  be  insignificant 
for  the  MOD  I  and  MOD  II.   The  small  effect  on  predicted  turbine  per- 
formance due  to  streamline  curvature  indicates  that  the  exact  value 
assumed  for  the  curvature  factor  K  will  not  greatly  affect  predicted 
performance.   The  results  of  this  analysis  also  indicate  that  Vavra's 
method  of  approximating  streamline  curvature  is  of  sufficient  accuracy 
for  methods  of  analysis  where  the  flow  equations  are  satisfied  at 
stations  between  blade  rows. 
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It  is  recommended  that  the  MOD  I  be  tested.  Results  of  that 
investigation  would  provide  additional  information  concerning  the 
general  applicability  and  accuracy  of  the  method  of  analysis  pro- 
posed in  this  thesis. 

Predicted  and  experimental  flowrates  for  the  MOD  II  were  in 
close  agreement.   Since  the  restriction  factor  £   is  a  significant 
factor  in  the  predicted  flowrates,  the  experimental  results  verify 
the  validity  of  the  theory  used  in  the  developement  of   £ 

Experimental  results  for  the  MOD  II  showed  that  the  predicted 
rotor  torque  and  turbine  efficiency  were  generally  2  to  3  per  cent 
too  high.   The  angles  and  velocities  predicted  for  the  flow  at  the 
stator  outlet  were  in  excellent  agreement  with  the  experimental  re- 
sults.  Therefore  the  high  values  predicted  for  rotor  torque  and  tur- 
bine efficiency  must  be  tied  to  the  rotor  solution.   The  measured 
values  of  outlet  angles  and  velocities  for  the  rotor  discharge  also 
indicate  that  the  predicted  rotor  solution  is  not  completely  correct. 

Exactly  what  parameters  are  in  error  in  the  predicted  rotor 
solution,  and  to  what  extent,  cannot  be  stated  without  more  experi- 
mental data.  More  traverses  should  be  taken  at  the  stator  and  rotor 
outlets  using  recently  calibrated  probes.  With  this  in format ion, loss 
coefficients  and  relative  discharge  angles  for  the  rotor  could  be 
determined  on  a  streamline  basis.   Based  on  the  experimental  results 
included  in  this  thesis, it  is  suspected  that  the  magnitudes  of  the 
predicted  discharge  angles  p2  are  too  large.   If  the  predicted 
turning  angles  of  the  flow  through  the  rotor  were  less,  the  predicted 
torque  and  efficiency  would  decrease.   Also,  the  separated  flow  at 
the  hub  and  tip  at  the  discharge  of  the  MOD  II  rotor  indicates  that 
predicted  loss  coefficients  should  be  higher  for  the  streamlines  at 
these  locations. 

To  accurately  predict  turbine  performance  at  off-design  as  well 
as  design  conditions, a  streamline  analysis  is  necessary.   The  experi- 
mental verification  obtained  for  the  proposed  method  indicates  that 
the  method  possesses  much  potential  and  that  additional  developement 
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is  warranted.   For  future  developement  work  on  this  type  of  analysis 
it  is  recommended  that  the  following  changes  be  considered: 

1.  Use  seven  streamlines  instead  of  five. 

2.  Calculate  discharge  angles   -   using  the  methods  of  this 
thesis  but  neglecting  blade  curvature  effects.   That  would  decrease 
the  angles  by  2  to  3  degrees  and  be  in  closer  agreement  with  experi- 
mental results. 

3.  Introduce  a  multiplying  factor  for  the  rotor  loss  coeffici- 
ents which  would  more  accurately  account  for  the  separated  flow  at 
the  hub  and  tip  but  would  not  greatly  change  the  average  rotor  loss 
coefficient.   For  example,  rotor  loss  coefficients  would  first  be 
calculated  using  the  methods  of  this  thesis. and  then  the  loss  co- 
efficients for  streamlines  1,  2,  6,  and  7  would  be  multiplied  by, 
say,  1.4  and  those  for  streamlines  3,  4,  and  5  would  be  multiplied 
by  0.6  or  0.7. 

4.  Calculate  rotor  loss  coefficients  for  an       range  of 
-3.0  to  2.0  instead  of  -2.0  to  1.6.   The  increased   —  range  would 
improve  performance  prediction  for  turbines  with  untwisted  rotor 
blades  at  off-design  conditions. 


49 


XI- 


50 


•w- 


-^ 


(-) 


VELOCITY      DIAGRAMS 
FIG.  3 


52 
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FIG.  4 
STATOR     a    ROTOR     BLADE     PROFILES      (MOD  I). 
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NOTE. 


BLADES    ARE    OF 
CONSTANT    SECTION. 


RADIAL     PLANE 
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SCALE:  2. 25  X 


FIG.  5 
STATOR    8    ROTOR     BLADE     PROFILES    (MODH). 
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APPENDIX  A 
DEVELOPEMENT  OF  EQUATIONS 

1.   Equation  of  Motion 

The  equation  of  motion  for  relative  flow  is 

\  I  R  =  \\.  x  (V  >  V-  +2  CD)  +  TVS 


In  cylindrical   coordinates: 


VHf 


I 


aHR 


■ 


(2) 


(87) 


W  *   2  UJ    =  tr  '  ■■    *>  VVg  -  Le  (2  a;  Wr) 


l      £S 


as 


o:  \ 


Equating  Eq.    87   to  the  components   in  Eqs.    88-90, 


.  9  ^  2       I 


(88) 

(89) 
(90) 

(91) 
(92) 

(93) 


11 


Eckert,    o£.    cit_. ,    pp.    149-155. 
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For  assumed  axisymmetric  flow,  — — —  «  0, 

Wa  Ur\Nu)  _  Wr    VrWu)  _  «       ^  (94) 


la.  0  = 


r    ^     r   ^r-a^wr 


£2  ^2 


7,  SHR    *A,    *Wr     >A/     3Wft    ,    Wy 


(95) 


|k.=  ^^^.WA||r+WA^  +  2a;Wu+T|^  (96) 


Equation  94  can  be  written 

Replacing    —   :     in  Eq.  95  by  its  equivalent  in  Eq.  97 


—   '  —  -J7--Zujr*W  <97> 


Hk-W     i-^-l      i^       >/Vu  /    WrSfrrWJA    ^^WtfWrtT^.  (98) 


or 


j^.WjWr  Vv»M  ^w  jWf,^^^  ^— ^  .T^,  (99) 


Multiplying  Eq.  96  by  W  and  Eq.  99  by  W  , 


^-^•iiLi._.^.v/v^r,VVfiWr^.  .,,..,  .WrT^  (100) 


•-if     r-    o^        t>2"       c  i 


Adding  Eqs.  100  and  101, 


rWr|^n) 


(101) 


(102) 


For  adiabatic  flow  IL  is  constant  along  a  streamline,  therefore, 

V\Mt  -7Hf  =c  (io3) 

For  axisymmetric  flow, 

(104) 

(105) 

V^-^(-^^)=T(wr||^l^         (106) 

^:   _  lli.   il  (107) 


From  Eq.    104, 

^HB   _     VW     3>_Hr 
■7r~  '   "   Wfl       -; 
Substituting  Eq.    105   into  Eq.    102, 
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Using  Eqs.  105  and  107,  Eq.  99  becomes 

■'-   -'*.-  L_l_        \\ir^-\Nr^-;uj^--^-T^  (108) 

Multiplying  Eq.  108  by  f  -  V-  A  )  >  gives 

^Hr   „  VV_  ^rM  _  yW  SJQ&  ..  ,■  .  ^i>  ,  ■  ..  I'V,  - "  ?7         (109) 

Eq.  109  is  identical  with  Eq.  96, therefore  Eq.  109  is  the  equation 
that  must  be  solved.  With  \'--'fi  -r-r  =   ~  '   [jEq.  109  becomes 

Rewriting  Eq.  110  for  the  rotor  discharge  and  using  the  substitutions 

M  sH    ^     and    T  =Jk.  W£.  ) 


(110) 


Cp  jte   ° 


(111) 


(112) 


The  last  term  of  Eq.  Ill  can  be  expressed  by 

With  /'—.--,    ^   and/|+TAM  h-      -\\Vq.  112  becomes 

Using  Eq.  113,  and  after  rearranging,  Eq.  Ill  becomes 


(114) 

cP  [       J  ^ 

Noting  that 

5  ( 

— --  -    —    '    -  (115) 

Eq.  114  is 


(116) 
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Th 


,  r„ 


is  equation  is  made  non-dimensional  by  multiplying  by  (  — -£--   )  > 


V/a! 


where  the  subscript  m  refers  to  the  mean  streamline.   In  the  sub- 
sequent derivations  the  subscript  2  will  be  omitted.  Then 


.4.  w.O 

A^_ 

^t 

q  r 

a/r 

^ 

Now  let 

x= 

r 

Equation 

-  -JL 
117 

then 

becomes 

W/£   WA;  Ur 


O  (in) 


(118) 


(119) 
(120) 


s^_i  .?l\  ^    XL  is^  +  2  y  ™ ' :  z(mtm/2) 


:>x 


+  + 


£>X 


9-"Mr-Y2T/Mz/S 


\Nt 


The  fourth  term  in  Eq.  121  can  be  expanded  and  rearranged  as 


O   (121) 


^T! 


r(XYTON,.-;);  =  ^Y^ 


■%2YTflNV|f  ^TAN/*  (122) 


^  A 


with 


and 


•/ r).-r        i         o/£ 


^y  C:./^     £>x 


2Y„„Vii.™„V4p 
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Equation   122   then   is 

aft*) 

^X 


Y2,   *Wr     Ya   ^>s* 


2T/wl  o/f  ,  2Y 


Multiplying  by  '  c-J-fi/s(  2 )  >  and  with 


TAM 


P 


+  TAN) 


il 

^x 

1 


o-/? 


Y: 


^^^^-*^^!f*f-^ 


(122a) 


w^Y 


The  terms     coe?/^ 


W^Y*   W^Y^X   W^Y 


and 


Wa  *a 

streamline  curvature  and  can  be  approximated  by 


represent  the  effects  due  to 


(122) 


(124) 


COS2^  = 


(125) 


L  +(^-) 
If  the  streamline  curvature  in  meridional  planes  is  zero,  the  terms 

represented  by  Eqs.  124  and  125  will  take  on  values  of  zero  and  one, 
respectively.  With  ir  positive  as  shown  in  Fig.  2,the  plus  sign 
is  used  for  K  at  station  2  and  the  minus  sign  at  station  1.   The 
streamline  slope  as  represented  by  Eq.  125  is  the  same  at  stations 
1  and  2  since  the  streamline  pattern  is  assumed  to  repeat  itself  after 
station  2. 

Using  Eqs.  124  and  125,  Eq.  114  then  becomes 


-—-cos, 


2Kf. 


m&) 


Z TAM/? 


t-T^ 


WAvnY  Wa^       ^/2  ^  1  IMLY*  ~s,^JclX       (126) 
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For  the  enthalpy  terms  to  be  in  terms  of  BTU/lb  ,they  must  be 

m 

divided  by  2gJ  to  keep  the  equation  dimensionless.  With 
C  -  2g  J,  Eq.  126  is 


W*„Y 


Noting  the  corresponding  terms  for  the  absolute  flow  of  the 
stator,  and  with  U=0,  the  equation  for  station  1  is 


dM3. 


dX 


CoS: 


4KSK-^)tf]— *4 


^ 

X 


2    <lklVl    *    C  cos;c*. 
-—  SIM  O.  +  - 


YlVA*     dX 


2.  Flow  Function  J 


Flowrate  in  lb  /sec  for  the  expansion  process  shown  in 
Fig.  44  can  be  written 


W  = 


r 


AV 


(i27y 


(128) 


(129) 


(PLENUM) 


A>p>T>V>/° 


W 


Expansion  From  Plenum 
Fig.  44 
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Assuming  an  isentropic  expansion, 


rt  \n< 


P  \T 


(130) 


and, 


VT.Tj,  [.-(f)*] 


(131) 


then  using  the  thermodynamic  relations 


V-  -T-T 


-p    J   tf-l 


(132) 


the  discharge  velocity  can  be  expressed  by 

L  L  ^ 

The  density   /^   is  found  with  the  isentropic  assumption  by 

Pt     P 


7^  = 

or 

/  P  \  7      fP\1  f? 

Using  Eqs.  133  and  135,  Eq.  129  becomes 


Hn^e 


'   >9^Ttj.(-(^-]( 


This  can  be  rearranged  to  give 


&[»)*- (*)11 


(133) 


(134) 


(135) 


(136) 


(137) 


A    Pr 


tftl    -,\l 


hmwr-Hrrt-i 


(138) 
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3.   Restriction  Factor   |   2 

The  expansion  process  through  a  blade  row  la  shown  in  Fig.  45 
(a).   The  letter  u  represents  the  velocity  at  some  point  within  the 
boundary  layer  and  V  .  represents  the  velocity  for  an  isentropic 
expansion.  The  various  dimensions,  angles,  and  velocities  for  the 
throat  location  are  shown  in  Fig.  45(b). 


/  / 


u 

-3^P 


T,  * 


23JCP 


'th 


(A) 


Conditions  at  Exit  of  Blade  Row 
Fig.  45 


where: 


The  mass  flowrate  per  unit  blade  height  can  be  written 


/^h  «  density  at  P.  for  isentropic  expansion  (slugs/ft  ) 
1°   -  density  corresponding  to  the  velocity  u  (slugs/ft  ) 

<Xd*  ■  blade  angle 
S  -  boundary  layer  thickness 


vavra,  M.  H.,  Problems  of  Fluid  Mechanics  in  Radial  Turbo- 
machines  (Rhode-Saint-Genese,  Belgium:  Von  Karman  Institutdcfor 
Fluid  Dynamics,  1965)  VKI  Course  Note  55b,  pp.  G46-50. 


(139) 


96 


For  Mach  numbers  of  0.8  or  less,  the  discharge  angle  can  be 
calculated  quite  accurately  by 

"1 


where: 


a  ■  minimum  throat  width 
t  -  blade  thickness 


then  with  Y,   -  —     and  using  Eq.  140  in  Eq.  139, 


which  can  be  expressed  as 


-=A^4«-^('-/Mv?)] 

Assuming  constant  pressure  at  the  throat, 

/   RT 

and 

where 

fth  '   T 

-y  -      ,.2   "     /  x/.?  \  /  *,  V*       w  /  u  \ 

In  Eq.  144  the  term  X  is 


<th 


i" 


(140) 


kV=A  ^IZM  cos.;   -2S(aS^  J  ^A^V/?*  ^  (     (     } 


(142) 


(143) 


(144) 


The  boundary  layer  is  assumed  to  be  turbulent  and  the  profile  can 
therefore  be  expressed  by 

t 


(146) 
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Then  with  the  displacement  thickness  of  the  boundary  layer  defined  as 


s*ss  i-o-x 


0 


fcf 


(147) 


.  142  can  be  written 


£«/foV<**[,~  "5"] 

The  loss  coefficient  to  the  throat  is 


M--fc 


(148) 


(149) 


th 


The  loss  coefficient  can  also  be  expressed  in  terms  of  the  kinetic 
energy  lost, by 

where  E  represents  energy  rate  at  the  discharge  due  to  the  average 
velocity  V,,  or 

#■        / 

E  -/%  Vf h  [*-  £<]  -^  +  £/("/  -r  ) d/  (I") 


Equation  151  can  be  rewritten  to  give 


E=Al«i-£i 


Then  using  the  energy  thickness  which  is  defined  as 
the  loss  coefficient  can  be  expressed  by 

>«■     -       -i-   ... 


^a 


3T- 1 


£^ 


(152) 


(153) 


(154) 
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The  restriction  factor  t  represents  that  part  of  the  throat 
opening  in  which  would  occur  the  uniform  theoretical  velocity, 
therefore 

l-l-^-il  (155) 


Now  defining  an  energy  parameter  H***  by 
1 


H**  =  -2—  (156) 


Eq.  155  becomes 


H***-l 


5  "  H***-i^ 

The  loss  coefficient  $      in  Eq.  157  accounts  for  the  losses  that 
occur  from  the  inlet  to  the  throat  of  the  blade  channel.   No  raeans 
exists  for  predicting  this  loss  coefficient.  Half  the  total  loss 
coefficient  for  the  blade  row  -y-  will  be  used  to  represent  these 
losses.   Eq.  157  then  becomes 


(157) 


5     H     -\+"2i 

*    13 
4.  Method  of  Evaluating  H*** 

I  use  of  the  binomial  theorem,  the  denominator  of  the  integral 

part  of  Eqs,  147  and  153  can  be  expanded ^yielding 

3      4™.  4-      &m 


(158) 


(J -Kef      J  =  UXe.^        *Xe    *         +Xe7       +  Xe    tf 
j  integral  of  Eq.    153  is  then 
(X    3»i.    .         s^         .    2.      7#n     .,  3     9»*     ^  4      ilw*.  "1    j 


Then, integrating  and  evaluating  gives 

p!2   j,.-i_,J^  +  i^.  +  *i!.,**i 

Equation  153  can  now  be  expressed  as 

£  '       ^5^+1         S*M.+-\        7^*1        9wa+)      \l«*+l  /      3m+|        S~+>        7»~H 


13Eckert,   op_.    cit . ,   pp.    159-160. 
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r** 

s 

also, 

k.  e       /n^4|      Swt+I       7^+1         SW+I       ll^  +  l 

/x  _,V_L_  +  _L_  +  Jk_+Jkl  ,U^t..\ 

In  a  similar  manner,  the  expression  for  Eq.  147  is 

The  energy  parameter  can  now  be  expressed  by 

J-  *.  J_  +lf_  +  *i!  ,  If!    *** 
u***_  X&-)   5>*  +  i   sj*t»j   7**t»)   e^i   n»*  +  i 

Xe     I         M+»  fm+l         fm-tl        7*i+l        e^-n 

The  number  of  terms  in  the  numerator  and  denominator  of  Eq.  159 
is  considered  sufficient  to  give  good  convergence  for  H***. 

Although  the  exponent  m  for  the  turbulent  boundary  layer  is 
dependent  on  Reynolds  number,  it  is  taken  to  be  a  constant  for 
this  analysis. 


(159) 
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APPENDIX  B 
COMPUTATION  OF  OUTLET  ANGLES  AND  LOSS  COEFFICIENTS 

1.  Outlet  Angles 

The  absolute  discharge  angles  for  the  stator  and  the  relative 
discharge  angles  for  the  rotor  are  computed  by  using  the  same 
methods.   For  these  calculations,  the  assumption  is  made  that  the 
outlet  angles  are  not  influenced  by  the  flow  incidence  angles. 
The  factors  which  do  affect  the  discharge  angles  are: 

1.  Blade  geometry;  this  results  in  the  angles  being  a  func- 
tion of  the  radius  since  spacing  and  possibly  profiles  change  from 
hub  to  tip. 

2.  Radial  tip  clearance;  the  effects  of  tip  clearance  are 
assumed  to  influence  the  rotor  flow  from  the  mean  streamline  out 
to  the  tip  with  the  largest  effect  being  near  the  tip, 

3.  Exit  Mach  number;  in  accordance  with  the  experimental 
results  surveyed  by  Ainley,  the  outlet  angles  are  a  function  of 
Mach  number.  Values  of  the  angles  are  calculated  for  Mach  numbers 
M  of  0.5  and  1.0.  Then  a  smooth  curve  of  outlet  angle  versus  M 

is  drawn  between  these  points  with  an  inflection  point  at  M-0.75. 
Below  h~0.5,  the  flow  angles  are  assumed  to  be  equal  to  the  value 
computed  for  M»0.5. 

Vavra's  formula  is  used  for  the  first  approximation  of  the 
stator  outlet  angles  for  M  ^0.5. 


where 


The  effects  of  blade  curvature  are  next  taken  into  account 
by  using  the  method  given  by  Ainley,^*" 

(X=oC*+4(-§-)  (162) 


14 
Ainley,  oj>.  cit. ,  pp.  3-4, 
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In  Eq.  162,  e  is  the  mean  radiu:  of  curvature  of  the  upper  surface 
of  the  blade  profile  between  the  passage  throat  and  the  trailing 
edge.  This  quantity  is  approximated  by 


e  = 


fib 


where  j  and  b  are  shown  in  Fig.  46. 


(163) 


Blade  Geometry 
Fig.  46 

The  discharge  angle  for  M=1.0  is  found  in  a  manner  similar 
to  that  described  by  Ainley, 


(X,    v  =  COS 

(v uo) 


fe) 


(164) 


The  change  in  0(   due  to  the  increase  in  local  Mach  number  from 
0.5  to  1.0  is  assumed  to  be  constant  along  the  blade  height  and 
equal  to  that  computed  for  the  mean  radius. 


$oCT  =  SocH=o(1 


oc, 


Relative  discharge  angles  for  the  hub  and  mean  radii  of  the 
rotor  are  found  by  the  same  method  as  used  for  the  stator.  The 
effects  of  tip  clearance  are  accounted  for  in  the  computation  of 
A^  for  the  tip  position.  With  Ainley1  s  formula,  the  flow 

angle 


(165) 


A 


for  the  tip  is 


a  -  ™4  -  ^feg)><  *  xfflteg- w-3 


(166) 
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The  new  terms  in  Eq.  166  are: 

h  -  blade  height 

/£a-  discharge  angle  before  tit  clearance  is  accounted  for 

X  -  factor  which  depends  on  the  type  of  shroud 

Ainley  suggests  that  X  be  set  equal  to  1.35  if  no  shroud  is 
arranged  at  the  rotor  tips.   However,  Ainley1 s  value  has  been 
determined  with  a  mean  streamline  analysis.   To  localize  the 
effects  of  tip  clearance,  X  has  been  changed  to  5.4,  a  value  which 
gives  about  the  same  overall  tip  clearance  effect  as  if  the  value 
of  X=1.35  were  used  for  the  mean  streamline.   Fig.  47  shows  that 
the  areas  under  the  two  lines  representing  X  as  a  function  of  radius 
are  approximately  equal.   The  slope  of  the  solid  line  representing 
X  for  this  method  is  obtained  because  of  the  assumed  linear  varia- 
tion of  outlet  angle  with  radius  between  the  mean  radius  and  the  tip, 


1.35 


Tip  Clearance  Factor  X  as 
Function  of  Radius 

Fig.  47 

Sample  calculations  for  stator  outlet  angles  are  shown  in 
Table  I.  A  sample  of  rotor  outlet  angle  computations  can  be  seen 
in  Table  II. 
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TABLE  I 

SAMPLE  CALCULATIONS  FOR  STATOR  OUTLET  ANGLES 

MOD  I  turbine 

M^0.5 
radius  -  in. 
spacing  (s)  -  in. 
throat  opening  (a)  -  in. 
trailing  edge  thickness  (t)  -  in. 
a 

(Sioof3  =A 

Kt=l-(0.0027)    (A)    (|)  0.9849     0.9891       0.9932 

a 
sK* 


3.597 

4.125 

4.950 

1.738 

1.994 

2.392 

0.455 

0.590 
-  0.044 

0.804 

0.2619 

0.2960 

0.3360 

21.4 

13.6 

7.48 

*<*> 


-1,  a  x 

cos    (iic^ 


.* 


0.2658  0.2991 

0 

.3384 

74. 6e   72.6* 

6 

70.2 

0.703   0.9106 

1.37 

75.3     73.5 

71.6 

M=1.0 

<X     »  of,  +  4(-) 


Use  4.125  in.  from  centerline  as  representative  radius. 

Projected  trailing  edge  thickness  (t*)-in.  0.13 

°W  -  cos"1  (^jf)  -  deg.  71.5° 
$<*   .  73.5  -  71.5  =  2.0 

°^  73.3    71.5     69.6 
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SAMPLE  CALCULA'i.  ,<NS  FOR  ROTOR  OUTLET  ANGLES 


MOD  II  turbine 


tip  clearance  =0.033  in, 


0.5 


radius  -  in. 

3.3235 

4.197 

4.918 

spacing  (s)  -  in. 

1.595 

1.4643 

1.7158 

throat  opening  (a)  -  in. 

0.382 

0.560 

0.721 

trailing  edge  thickness 

(t)  -  in. 

0.048 

a 
s 

0.3295 

0.3824 

0.4202 

!)Goo) 

4.1397 

3.278 

2.7975 

&-.-T.A 

109.0 

50.3 

29.8 

lk*|-<0.oo*7)(A)(f.) 

0.9030 

0.9481 

0.9662 

a 

0.3649 

0.4033 

0.4349 

**-«'(^ 

-68.5° 

-66.2° 

-64.2 

*(*) 

1.75 

2.37 

3.04 

/s-/2*-*m 

■   -70.3° 

-68. 6° 

-67.2 

k    -  in. 

0.033 

*  /h 

0.0207 

X-(4)  (1.35)=type  shroud 

factor 

5.4 

Cos^0  ,  ^=50° 

0.86603 

Cos^',^^-^7.20 

0.38268 

«&)£*•« 

0.25297 

tan  ^ 

-2.4142 

(l-B)tan/4' 

-1.8035 

tan  X> 

0.57735 

(l-B)tan^\BT/IM^  : 

-c 

-1.6574 

/^-tarT^C) 

-58. 94 
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TABLE  II  (continued) 
SAMPLE  CALCULATIONS  FOR  ROTOR  ANGLE 

M2«1.0 

Use  radius  «*  4.197   in.   as  representative    radius. 
Projected  trailing  edge  thickness   (t*)-in.  0,108 

&^C°**(^T*)  -65.6 

S/f  —  68.6+65.6=  -3.0 


radius   -   in.  3.3235       4.197  4.918 

■62.6° 
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2.   Loss  Coefficients 

Loss  coefficients  were  predicted  using  Ainley's  methods. 
The  techniques  used  are  completely  described  in  Refs.  1  and  2. 
The  assumptions,  definitions,  and  equations  which  are  necessary 
for  a  basic  understanding  of  the  method  are  described  below. 

It  is  assumed  that  loss  coefficients  are  independent  of  Mach 
number  and,  for  a  given  glade  geometry,  are  a  function  of  the  flow 
incidence  on  the  blade.   Losses  are  divided  into  profile  losses, 
secondary  losses,  and  losses  due  to  tip  clearance.  Mixing  losses 
are  not  accounted  for  separately  but  are  included  in  profile  and 
secondary  losses. 

A  parameter  Y  is  defined  by 

Y  =  Loss  in  total  pressure 

Total  pressure  at  outlet  -  outlet  static  pressure         (167) 

For  incompressible  flow,  which  is  assumed  for  these  calculations, 
the  loss  coefficient  is 

<*/  =  ^  t  (168) 

J      \+£1j 
The  subscript  y  on  Y  in  Eq.  168  indicates  the  type  loss  repre- 
sented; i.e.,  Y  would  be  the  profile  loss. 

For  zero  incidence,  the  Y  for  profile  losses  is 

where: 

^  -   blade  discharge  angle  not  accounting  for  tip  clearance 

t  -  maximum  blade  thickness 
m 

c  -  blade  chord 

Equation  169  is  applicable  to  both  the  stator  and  rotor.   However, 

for  stator  use  the  blade  inlet  angle  j£q     is  zero.  The  quantities 

io/j,      fl/vand  \p/a    Aare  taken  from  Fig.  4  in  Ref.  2  where  these 

quantities  are  plotted  as  a  function  of  blade  solidity.   The  maximum 

value  of  t  ,   is  set  equal  to  0.25  for  thick  blades, 
m/c        M 

To  determine  profile  losses  for  incidences  other  than  zero,  Fig. 

6  in  Ref.  2  shows   3 p    versus  4—  where  i  is  defined  as  that  in- 

Yp«--o)        1%  s 

cidence  where  the  losses  are  twice  the  minimum  losses. 
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Secondary  and  tip  clearance  losses  are  computed  by 


where: 


/^L-  the  mean  of  the  gas  inlet  and  exit  angles  of  the  blade  row 
6  -  factor  representing  type  shroud,  similar  to  X 
^-a^K^-TAM^Cos^-  fit) 


a-^(^)7(.+  ^)] 


(171) 


in  Eq.  171: 

A,  ■  (annulus  area  at  discharge)  (cos  /&^      ) 

A.  «  (annulus  area  at  inlet)  (cos  ^     ) 

I.D.  and  O.D.  -  inner  and  outer  diameters,  respectively,  at  the 
blade  exit 

Values  of  f)       are  obtained  from  Fig.  8  in  Ref.  2.   Ainley  suggests 
a  value  of  0.5  for  B  for  the  type  shroud  used.  With  the  same  rea- 
soning as  previously  used  for  the  effect  of  tip  clearance  on  dis- 
charge angles,  B  is  set  equal  to  2.0  for  the  losses  at  the  tip  and 
zero  for  the  mean  radius  and  hub  locations. 

Sample  calculations  for  the  prediction  of  stator  loss  coeffi- 
cients are  shown  in  Table  III.   Sample  calculations  for  prediction 
of  stalling  incidence  and  rotor  loss  coefficients  are  shown  in 
Table  IV. 
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I  BLE  III 
SAMPLE  CALCULATION  FOR  STATOR  LOSS  COEFFICIENT 


MOD  I  turbine 

radius  =  3.597  in. 

pitch  (s)  -  1.738  in. 

All  figures  referred  to  are  in  Ref.  2 


inlet  angle  (  <X0  )  ■  0 
outlet  angle  (  <X,  )  -  75.3* 
chord  (c)  =  2.691  in. 


\=o^)\^[$) 


Y  -  Y 


VA-Voc, 


<XQ-o 


K'-*)   =  °'^7 


Prrcrt   __  1730 
CHOKZ>'  2.631 


=  o,64S9 


"  - '  [Wo * H)]  ■  *  [^  7"')l/c>  ♦  &yl 

h=  { (0.038L)  h-  O.OOS4,  (f^8) 


O,  D,  -  9. 538  4/1. 


c. 


~  =  £  [tAn>  0to  -TA/0  OL^j  cos  *  ^  -  2  [-  (TAN)  -75.3°)]  Co*  £?.  1° 


%■  =  3.543^ 


Y«  = 


ro.oos6ir3.s458ia[>^7S-s8~U  o.o 


45  1 


Y*  +  Yp  =  o.o9a9 

CORRECTION     FACTOR      FoR     Tfcfi\UW6    ED6E    THICKMESS    ,    F=1.03 
(j.os)(o.  0929s)  =  00949 


Q-Q949 
\-t-0. 094-9 


=    0.06<o(o 
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TABLE  IV 

SAMPLE  CALCULATIONS  FOR  STALLING  INCIDENCE  AND 

ROTOR  LOSS  COEFFICIENTS 


MOD  II  turbine 
radius  ■  4.918 


t    =  0.728  in. 
max 


s  -  1.7158  in. 
k  -  0.033  in. 
c  «■  1.967  in. 


A   -  30° 
A'-  -67.2° 

blade  height(h)  -  1.595  in. 


Figures  referred  to  are  in  Ref.  2 

lA^ 


I-  Q-7Z& 


C 


s 

A 


■A')"V-°) 


%_\"%' 


7"  <«■  »")  ,  -  A  .-0.444  ;   fc^*.  o.l393 
-   '-7'gg.      .-.  .,.;  Y  -o««*      y  „A=°'IS4       6% 4") 

">P;=     =?0.036S+^O.i9S3Xo.iS4-o.oiiS)?(|.2f")0'*4    -    0-O6,C>Z 
"X-  =  6.9SS    (Fig.7)       /$  --70.4*         —2-     =6,426 


(|-M*) 


'£ 


(i.Mi5"«-*''fc'»""7"  CFi9'7)   '    'ts5e.s°-7'=   2SJ 


VY,-^s(|)][y%]^£]     «...     8l.0iei.#1 

.         /-[/A    /    \2y         X  D  \1  r,D.  ootlET=  3-"52%5  L/l,        OX>.  ooTLET  s4,9lfl  lA. 

[>  v/\.+    O.P.  )J         *•*>,  IW»-ET  =  5,324S  Irt,         O».Itf|.ET-4.8S0€Ul. 

A  s  ir^4,9i82-g.5^3S2)    Cot  7o.i 
h=i(o,\48)  }  h-  0,0072.  (f^-fi) 

Y,  ^yk  =  [0,^2  +o.o)6ei][Vs/of [*3gy 
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TABLE  IV  (continued) 
SAMPLE  CALCULATIONS  FOR  ROTOR  LOSS  COEFFICIENT 


Y^ 

-1.5 

0.0 

1.5 

1-Cfc.Vi 

-38.25 

0 

38.25 

fit  =  i  +  /80 

-8.25 

30 

68.25 

T/M/f, 

-0.1450 

0.5774 

2.5065 

-1.2619 

-0.9008 

0.0638 

Aw=  TflN'7A) 

-51.6 

-42.0 

3.6 

*»*>£* 

0.6212 

0.7430 

0.9980 

tam/^-taK)/^  B 

2.2339 

2.9563 

4.8854 

%■*(.&)  ">»/~ 

2.7752 

4.3928 

9.7510 

NT 

7.7016 

19.2964 

95.0813 

'/co*V~ 

4.1726 

2.4385 

1.0061 

y^Y^(o.00i6^]^ 

0.2345 

0.3434 

0.6980 

VvP^  &*$ 

2.1 

1.0 

4.5 

v,  ■  lV%.e)]  V4 

0.1390 

0.0662 

0.2979 

YP*VYK 

0.3735 

0.4096 

0.9959 

(CORRECTION   F/KTOl5^(SgY)»C 

(%5)  0.3922 

0.4301 

1.0457 

3=  — 

0.2817 

0.3007 

0.5112 
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APPENDIX  C 
COMPUTER  PROGRAM 

A  computer  program  using  Fortran  IV  was  written  for  use  with 
the  performance  analysis.  The  program  was  used  to  predict  perfor- 
mance values  for  the  MOD  I  and  MOD  II  turbines.  However,  an 
attempt  was  made  to  keep  the  program  general.   If  the  methods  of 
this  thesis  are  utilized  for  expressing  turbine  characteristics, 
this  program  can  be  used  for  other  single  stage  axial  turbines. 

Input  for  the  program  consists  of  information  representing 
the  turbine  characteristics,  an  indicator  specifying  the  detail 
of  output  desired,  and  the  conditions  for  which  performance  values 
are  to  be  obtained. 

Twenty  input  cards  are  used  to  introduce  the  turbine 
characteristics.  The  information  contained  on  each  of  these  cards, 
with  required  dimensions,  is  listed  below. 

1.  Stator  mean  streamline  gas  outlet  angles  (ALF AM- radians) 
for  mean  streamline  exit  Mach  numbers  of  0.5,  0.7,  0.75,  0.8  and 
1.0. 

2.  Relative  rotor  mean  streamline  gas  outlet  angles  (BETAM- 
radians)  for  the  same  Mach  numbers  listed  in  (1.). 

3.  Radii  ahead  of  the  stator  (RC-in.)  for  the  five  stream- 
lines such  that  the  flow  area  is  divided  into  four  equal  parts. 

4.  Assumed  radii  of  the  five  streamlines  at  the  stator  exit 
(RS-in.). 

5.  Same  as  (4.),  only  for  the  rotor  exit  plane  (RR-in.). 

6.  The  predicted  stator  loss  coefficients  (ZETAS)  for  the 
hub,  mean  radius,  and  tip. 

7.  Ten  values  of  incidence  ratio  (RINC)  ranging  from  -2.0 
to  1.6  in  increments  of  0.4. 

8.  Ten  values  of  rotor  loss  coefficients  (ZETAR1)  for  the 
hub  corresponding  to  the  incidence  ratios  of  (7.). 

9.  Same  as  (8.),  only  for  the  mean  radius  (ZETAR3). 
10.   Same  as  (8.),  only  for  the  tip  (ZETAR5). 
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11.  Length  L  shown  in  Fig.  2  (CL-in.)  and  the  curvature 
factor  K  (CK) . 

12.  Differences  in  stator  gas  outlet  angles  (DALF-radians) 
from  the  angle  for  the  mean  streamline,  for  streamlines  1,  3,  and  5. 

13.  Same  as  (12.),  only  for  the  relative  rotor  outlet  angles 
(DBET) . 

14.  Ten  radii  for  the  stator  exit  plane,  equally  spaced  or 
approximately  equally  spaced,  ranging  from  the  radius  of  the  hub 
to  the  radius  of  the  tip  (Rl-in.). 

15.  Ten  values  of  stator  throat  opening  (Al-in.)  corresponding 
to  the  radii  of  (14.). 

16.  Same  as  (14.),  only  for  the  rotor  exit  plane  (R2-in.). 

17.  Ten  values  of  rotor  throat  opening  (A2-in.)  corresponding 
to  the  radii  of  (16.). 

18.  Inlet  blade  angles  for  the  rotor  (BETO-degrees)  for  the 
hub,  mean  radius,  and  tip. 

19.  Ten  stall  incidence  angles  (STALI-degrees)  corresponding 
to  the  radii  in  (14.). 

20.  The  radial  tip  clearance  of  the  rotor  (TIPC-in.),  the  number 
of  stator  blades  (ZS),  and  the  number  of  rotor  blades  (ZR) . 

Input  card  number  21  specifies  the  output  to  be  printed,  and  its 
use  will  be  described  later.   The  remaining  input  specifies  the  con- 
ditions for  which  performance  values  are  to  be  found, and  enters  the 
estimations  of  the  flow  Mach  numbers  used  for  the  first  approxima- 
tions in  the  iteration  process.   Input  card  number  22  specifies  the 
number  of  sets  (NSETS)  of  operating  conditions  for  which  solutions 
are  to  be  found.   A  card  containing  the  following  information  is 
used  for  each  point  specified  by  NSETS: 

a.  Estimated  Mach  number  of  the  flow  ahead  of  the  stator  (AMC). 

b.  Estimated  Mach  number  of  the  mean  streamline  flow  after  the 
stator  (AMS). 

c.  Total  inlet  pressure  (PTO-psi.). 

d.  Total  inlet  temperature  (TT0-°R) . 
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e.  Estimated  relative  Mach  number  of  the  mean  streamline 
flow  at  the  rotor  exit  (AMR). 

f.  The  rotor  speed  (RPM) . 

g.  The  ratio  of  total  inlet  to  static  discharge  pressure  (PR). 

There  are  eight  subroutines  in  addition  to  the  main  or  execu- 
tive part  of  the  program.  The  subroutines  and  their  main  functions 
are  listed  below: 

1.  Subroutine  PARAB  is  used  to  determine  the  coefficients  of 
the  parabolic  equations  used  to  approximate  curves. 

2.  Subroutine  LSQPOL  determines  by  the  method  of  least  squares 
the  coefficients  of  the  fourth  order  polynomials  used  to  approxi- 
mate curves. 

3.  Subroutine  CHAN  computes  the  flowrate  and  the  reference 
flowrate  from  conditions  ahead  of  the  stator. 

4.  Subroutine  STATOR  determines  the  axial  velocity  ratios  that 
satisfy  the  equation  of  motion  at  station  1. 

5.  Subroutine  FLOWREF  computes  the  reference  flowrate  at  the 
stator  exit  and  the  rotor  exit. and  adjusts  the  axial  velocity  of  the 
mean  streamline  flow  at  these  locations  to  satisfy  overall  continuity, 

6.  Subroutine  SLINE  checks  streamline  continuity  for  the  stator 
exit  and  the  rotor  exit, and  determines  new  streamline  radii  to  satisfy 
streamline  continuity. 

7.  Subroutine  R0T0R1  converts  the  absolute  flow  properties 
ahead  of  the  rotor  to  relative  flow  properties. 

8.  Subroutine  R0T0R2  determines  the  axial  velocity  ratios  at 
station  2  that  satisfy  the  equation  of  motion. 

The  main  program  computes  the  mass-flow-weighted  value  of  total 
inlet  to  static  discharge  pressure  ratio, and  adjusts  the  Mach  number 
of  the  flow  ahead  of  the  stator  to  obtain  the  pressure  ratio  speci- 
fied. The  overall  turbine  performance  values  of  tf   ,  M_,H.P.,r*,  and 
k.   are  also  computed  in  the  main  program. 

Computer  output  representing  the  solution  for  each  set  of  condi- 
tions is  printed  on  two  pages.   The  flow  properties  at  the  stator 
exit  plane  are  printed  on  one  page  .  '  A  sample  of  this  output  can 
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be  seen  in  Table  V.   The  second  page  contains  the  flow  properties 
for  the  rotor  as  well  as  the  overall  turbine  performance  values. 
A  sample  of  the  second  page  of  output  is  in  Table  VI.   Most  of  the 
output  quantities  are  self-explanatory.   However,  the  following 
symbols  may  not  be  readily  recognized  and  are  defined  as: 

SLINE  -  streamline 

ZETAP  -  one-half  the  total  loss  coefficient,  used  for  calcula- 
ting  | 

A  -  throat  opening  of  the  blade  channel 

W-FRAC  -  fraction  of  the  flowrate  w  between  the  hub  and  the 
indicated  streamline 

BO  -  rotor  blade  inlet  angle 

INCID  -  rotor  flow  incidence  angle 

W 
PSIR  -  relative  flow  velocity  coefficient    2 

W2th 
In  addition  to  the  main  output  there  are  additional  WRITE 

statements  in  the  program.   The  additional  output  allows  the  user 
to  follow  more  closely  the  intermediate  steps  of  the  solution  pro- 
cess.  This  output  is  als~v  helpful  when  program  changes  are  made 
which  require  debugging.   The  extra  output  is  not  desirable  for  pro- 
duction runs  because  of  tf-e  large  amounts  of  computer  running  time 
and  printout  that  result.   When  the  additional  output  is  wanted,  1 
is  read  into  the  computer  on  input  card  number  21  for  the  indicator 
IND.   If  0  is  read  into  the  computer  for  IND,  only  the  main  output 
representing  the  final  solution  will  be  printed. 

The  statements  of  the  program  are  listed  on  the  following  pages, 
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